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(57) Abstract: The present invention 
relates to a method for wavclcnglh 
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TPSF -based optical imaging system. 
This consists of identifying several 
chromophoccs in a highly rarbid 
medium and selecting optimized 
wavelengths whereby using these 
wavelengths optimizes the deduction of 
the chromophorc concentrations- Such 
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combined to deduce other properties of 
the turbid medium. 
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CHOICE OF WAVELENGTHS FOR MULTIWAVELENGTH 
OPTICAL IMAGING 

Field of the Invention 

The present invention relates to the field of optical imaging in which objects 
5 which diffuse light, such as some human body tissues, are imaged using signals 
resulting from the injection of light into the object and detection of the diffusion of 
the light in the object at a number of positions. More particularly, the present 
invention relates to the choice of wavelengths for multiwavelength optical 
imaging in order to provide enhanced information. 

in Background of the Invention 

Time-domain optical medical images show great promise as a technique for 
imaging breast tissue, as well as the brain and other body parts. The objective is 
to analyze the temporal point spread function (TPSF) of an injected pulse as it is 
diffused in the tissue, and the information extracted from the TPSF is used in 
15 constructing a medically useful image.- 

For example, one can extract time-gated attenuation information from the TPSF 
which improves the image spatial resolution over previous continuous wave 
methods. However, it is unclear whether such Improvements in image spatial 
resolution are adequate for diagnosing breast cancer based on morphology. 

20 An alternative approach is to use the TPSF to decouple the light attenuation into 
absorption and scattering components. This extra information.- which cannot be 
obtained from continuous wave methods, may be clinically useful. Moreover, 
one can obtain the tissue . absorption spectrum by performing time-domain 
measurements at multiple wavelengths. In tissue there are several molecules 

25 which absorb the light and are known as chromophores. Spectroscopic analysis 
of the tissue absorption spectrum peimits chromophore concentrations to be 
measured. Furthermore, combination of the chromophore concentrations can 
yield physiological information, as opposed to morphologic information, which 
could provide a more medically useful image. 

30 The problem is one of knowing which are the dominant chromophores to include 
in a tissue model and then choosing the "best" wavelengths to deduce their 
concentrations most accurately. 
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Summary of the invention 

It is an object of the invention to improve image quality in TPSF-based optical 
images by choosing an efficient combination of wavelengths and combining 
information from the combination of wavelengths. 

5 It is an object of the present invention to provide an objective method for 
choosing the wavelengths for a multiwavelength TPSF-based optical imaging 
approach. For a given set of chromophores, the best selection of the 
wavelengths is performed for the set as a whole as opposed to choosing the 
best wavelength for each chromophore individually. Furthermore, hardware 
10 constraints can be taken into consideration in order to optimize the selection of 
wavelengths for a given device. 

Brief description of the drawings 

Fig. 1 illustrates the absorption spectra used of oxy-Hb, deoxy-Hb, pure water 
and lipid; 

15 Fig. 2 illustrates the inverse of the condition number of the hemoglobin specific 
absorption matrix as a function of wavelength X1 being plotted for a system of a) 
two wavelengths where the other wavelength is fixed at X2 = 850 nm, b) three 
wavelengths where the other wavelengths are fixed X2 = 850 nm and \Z = 758 
nm and c) four wavelengths where the other wavelengths are fixed 12 = 850 nm, 

20 X3 = 758 nm and 7A - 800 nm; 

Fig. 3 illustrates the inverse of condition number C for the specific absorption 
spectra of oxy-Hb and deoxy-Hb as a function of X1 and X2. The plot is 
symmetric with respect to the diagonal. Regions of high values indicate 
combinations of wavelengths advantageous for spectroscopy; 

25 Fig. 4 illustrates the deviation of calculated saturation and true saturation 
(S(calc)-S(true)) for a model tissue containing 15 \iM [HbT], S(true) = 25%, 50% 
and 75% and a lipid concentration of 40%. Two wavelengths at 760 and 850 nm 
were used to fit [oxy-Hb] and [deoxy-Hb]. The sensitivity with respect to wrong 
assumptions of lipid and water concentrations are shown; 

30 Fig. 5A illustrates the inverse of condition number C for the specific absorption 
spectra of oxy-Hb and deoxy-Hb and lipid for a fixed wavelength X = 830 nm as 
a function of X1 and X2. The islands of high values indicate advantageous . 
wavelengths (scaling 0-0.01); 



01/16/2004 15:37 FAX 613 230 6706 



OGILVY RENAULT 



WO O3/0078OS> PCT/CA02/01081 

-3- 



Fig. 5B illustrates the inverse of condition number C for the specific absorption 
spectra of oxy-Hb and deoxy-Hb and lipid for a fixed wavelength X3 = 830 nm as 
a function of X1 and X2 (same as Fig. 5A but scaling 0-0.0005); 

Fig. 6A illustrates the inverse of condition number C for the specific absorption 
5 spectra of oxy-Hb and deoxy-Hb, lipid and water for two fixed wavelength at ^3 = 
760 nm and X4 - 830 nm as function of M and \2 (scaling 0-0.0015). Regions 
of high values are advantageous for spectroscopy; 

Fig. 6B Illustrates the inverse of condition number C for the specific absorption 
spectra of oxy-Hb and deoxy-Hb, lipid and water for two fixed wavelength at 7& = 
10 760 nm and X4 = 830 nm as function of M and \2. (same as Fig. 6A but scaling 
0-0.0005); 

Fig. 7 illustrates the estimation of deviations from true saturation values for a 
model tissue of [HbT] = 20 uM. S= 75%, a lipid concentration of 40% and true 
water concentration corresponding to 0-100% water. Three wavelengths at 760. 
15 780 and 850 nm were used for back calculation of S, shown here as a function 
of assumed water concentration; 

Fig. 8 illustrates the estimation of the influence of errors (noise) in ^ on the 
calculated Hb concentrations and saturation values. A model pa-spectrum based 
on 20 uM [HbT], S=50%, and a lipid and water concentration of 30% and 40% 
20 was assumed. Matrix inversion was performed for wavelengths 760, 790, 830 
and 850 nm. plotted is the change in calculated [oxy-Hb], [deoxy-HbJ and 
saturation value when the u« value at a single wavelength was changed by 
+0.0001 mm" 1 . This plot suggests that noise at 830 nm translates in the highest 
noise in saturation values; 

25 Fig. 9 illustrates the estimation of the recovery of saturation values based on 
different wavelength combinations. A model tissue of 20 pM [HbTJ, a true 
saturation of S=75%, lipid and water concentration of 40% were used, in the 
lower plot an offset of 0.0005 mm" 1 independent of wavelength was added to the 
model tissue pa-spectrum (no offset in the upper plot). Plotted are deviations of 

30 the saturation values due to matrix inversion and the true 75% value. The 
following wavelength combinations were used: 1) 760 nm and 850 nm, 2) 760, 
830 and 850 nm, 3) 760, 780, 830 and 850 nm, 4) 750-850nm, 5) 720-850 nm, 
6) 720-900nm; 

Fig. 10 illustrates the estimation of the recovery of saturation values based on 
35 different wavelength combinations. A model tissue of 20 uM [HbT], a true 
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saturation of S-50%, lipid and water concentration of 40% were used, in the 
lower plot an offset of 0.0005 mm" 1 independent of wavelength was added to the 
model tissue jia-spectrum (no offset in the upper plot). Plotted are deviations of 
the saturation values due to matrix inversion and the true 75% value. The 
5 following wavelength combinations were used: 1) 760 nm and 850 nm, 2) 760, 
830 and 850 nm, 3) 760, 780, 830 and 850 nm, 4) 750-850nm, 5) 720-850 nm, 
6) 720-900nm; 

Detailed Description of the Invention 

In accordance with the present invention, there is provided a method for 
1 0 selecting wavelengths for multiwavelength optical imaging. 

Tissue Chromophores 

The dominant near infrared chromophores contained in breast tissue are 
considered to be hemoglobin (Hb) in its oxygenated (oxy-Hb) and deoxygenated 
(daoxy-Hb) forms, water and lipids. Fig. 1 shows the absorption spectra of oxy- 

15 Hb (at 10jiM concentration), deoxy-Hb (at 10jiM concentration), pure water 
(100% concentration), lipid (absorption spectrum of olive oil has been used to 
estimate the absorption spectrum of fat). There are other interesting near 
infrared chromophores, such as glucose and cytochrome c oxidase, but their 
absorption contribution in the breast is considered negligible compared to the 

20 aforementioned chromophores. 

Physiological information 

Potentially useful physiological information about the breast tissue can be 
obtained from concentrations, Q, of the chromophores. The total hemoglobin 
concentration, [HbT], defined as [HbT] = [oxy-Hb] + {deoxy-Hb], is related to the 

25 local vascular density. Since cancer is commonly associated with an increase in 
vascularisation (angiogenesis), a measurement of [HbT] could be medically 
useful. The fraction of hemoglobin that binds to oxygen is known as the oxygen 
saturation, S, and defined as S=Ioxy-HbI/[HbT]. Increased metabolic activity 
increases oxygen demands which decreases the oxygen saturation. Since 

30 cancer is commonly associated with increased metabolic activity, a 
measurement of S could also be medically usefuL 

Wavelength Choice 

Historically as the biomedical optics field evolved the wavelengths were chosen 
for each chromophore individually by observing strong near infrared spectral 
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features for the given chromophore and using the closest hardware-available 
wavelength. Many researchers also used the isobestic wavelength of oxy-Hb 
and deaxy-Hb, the wavelength where their absorption per concentration are 
equal, since this wavelength is insensitive to the oxygenation state of the 
5 hemoglobin and can be related to the IHbT]. 

However, the question both posed and addressed here is that for a given set of 
chromophores what are the optimal wavelengths to use in order to deduce the 
concentration of each chromophore? It is interesting to note that the isobestic 
wavelength used by many researchers turns out not to be one of the 
1 0 wavelengths of choice. 

It Is an object of the present invention to provide an objective method for 
choosing the wavelengths for a multiwavelength TPSF-based optical imaging 
approach. For a given set of chromophores, the best selection of the 
wavelengths is performed for the set as a whole as opposed to choosing the 

15 best wavelength for each chromophore individually. Moreover, h is also possible 
to investigate scenarios such as the influence on determining chromophore 
concentrations under certain assumptions about the concentration(s) of other 
chromophore(s) in the set. Furthermore, hardware constraints can also be taken 
into consideration in order to optimize the selection of wavelengths for a given 

20 device. Fortunately, the recent advent of turn-key, pulsed, tunable near infrared 
wavelength lasers has permitted more viable availability of near infrared 
wavelengths. 

Experimental brute force approach 

One possible approach to optimize the choice of wavelengths for a given set of 
25 chromophores is to conduct a brute force experimental study. This would consist 
of performing numerous experiments where different combinations of 
wavelengths are evaluated for the given set of chromophores at known 
concentrations until the optimum combination for deducing their concentrations 
is found. Obviously, this approach' is likely to be highly time-consuming and it is 
30 not always trivial to provide a set of chromophores at known concentrations, 
particularly in the case of In vivo breast tissue. 

Matrix Inversi n sensitivity approach 

An alternative approach which avoids the numerous experiments of the 
experimental brute force approach is a matrix inversion sensitivity approach. 
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The equation which needs to be solved can be written for each wavelength as: 

J 



5 MM = T. m M> C ' 

t 

where ^ is the measured absorption coefficient, m a is the specific absorption 
coefficient of the different chromophobes and q is the corresponding 
concentration. 

1 0 This is written is matrix form as: 

lia = M-c 

where printing in bold indicates a matrix or vector, is a vector with a number of 
rows corresponding to the number of wavelengths (nO- c is a vector with the 
number of rows corresponding to the number of chromophores (n c ). Nl is a 
1 5 rectangular matrix of size n x x n<> 

If n x = n c the system can be solved by matrix inversion c=M" 1 p Q and if r\\ > n c the 
system is overdetermined and can be solved by the pseudo-inverse M* = (M T 
M)" 1 Wl T where M T is the transposed matrix of M. 

c = (M>a 

20 The pseudoinverse M* is an n x x n c array which is unique. If M is square (i.e. not 
overdetermined) the M* = WI" 1 . For given (i.e. chosen) wavelengths the 
pseudoinverse M + can be precalculated once and the matrix inversion 
corresponds to a simple matrix multiplication. This is the basis for the calculation 
of chromophore concentration. 

25 One means to quantify the expected sensitivity of a matrix inversion of a matrix 
M with respect to small errors in the data is the condition number C which is 
defined as: 

C = norm (M) • norm (WI" 1 ) 

C gives an indication of the accuracy of the r suits and is an estimate of the 
30 cross-talk between the different channels (i.e. chromophores concentrations)- 
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Values of C near 1 1ndicate a well-conditioned matrix. large values indicate an ill- 
oondrtioned matrix. The condition number is closely related to singular value 
decomposition (SVD) as rt is the ratio of the largest and the smallest singular 
value of a matrix. 

5 The matrix M for oxy-Hb and deoxy-HB at X = 780 and 770 nm is 

_ 0.3871 0.1465 X - 760wn 
"0.3280 0.1625^ = 770™ 

A matrix inversion is possible as the rank (M) = 2 F however the absorption at the 
two wavelengths is 'similar'. The condition number is C = 20.49. Choosing the 
wavelengths to be X = 760 and 850 nm gives the matrix 
HA w 0.3871 0.1465 

0.1729 0.2645 

Inspection by eye already shows that the absorption is very 'different'. This is 
confirmed by the condition number: 03.206. In what follows below the inverse 
of the condition number is plotted and analyzed. It has value between 0 and 1 „ 
1/C close to 1 means •orthogonal 1 spectra and low sensitivity to cross-talk. Small 
15 values of 1/C mean an ill-conditioned matrix. To find the best wavelengths, 1/C 
is calculated as a function of a wavelength. The wavelengths that give the 
highest values of 1/C are the best for a calculation of chromophore 
concentrations and the subsequent physiological information such as oxygen 
saturation, S. 

20 Model absorption spectra were generated with the absorption spectra of Fig. 1 
based on estimations of [HbT], S, lipid and water concentration. Matrix inversion 
based on different sets of wavelengths were performed to recover these 
parameters. These parameters were compared wtth the true ones for the 
different wavelengths and the sensitivity to noise or measurement offsets 

25 considered. 

Assuming that we. fit for the hemoglobin concentrations only and assuming 
certain values for water and lipid concentration, for a x-wavelengths matrix 
inversion, the best combination of wavelengths to give a well-conditioned matrix, 
the sensitivity of calculated values of oxy-Hb and deoxy-Hb concentration and 
30 oxygen saturation for variations of lipid or water concentration and sensitivity of 
S to measurement noise have been determined. 
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In Fig. 2 the inverse of the condition number is shown for matrices of oxy-Hb and 
deoxy-Hb specific absorption coefficients for 2, 3 and 4 wavelengths. In each 
case one wavelength (Xi) was varied between 650 and 950 nm while the 
remaining wavelengths were fixed 7* = 650 nm (2-wavelength system), %z = 850 

5 nm and ^ = 758 nm (3-wavelength system), and X* = 850 nm, X3 = 758 nm and 
JLj = 800 nm (4-wavelength system). Fig. 2 indicates that the selection of two 
wavelength at h\ = 850 nm and X 2 = 700 nm gives the highest values of 1/C and 
when the wavelength range is restricted via hardware constraints to > 750 nm, a 
system that includes the peak wavelength of deoxy-Hb close to 760 nm is 

10 advantageous. It does not matter whether two or more wavelengths are used. 
This somewhat counterintuitive result is valid only, without measurement noise 
and noise in the background absorption. 

Fig. 3 further highlights this finding for a two-wavelength matrix inversion. In this 
figure 1/C is plotted as a function of both at X1 and X z in the range 650-950 nm, 
15 The plot is symmetric with respect to the diagonal. Regions of high 1/C-values 
can be chosen and the corresponding 'good 1 wavelengths can be read off the 
axis. It is apparent that (with the restriction to > 750 nm) the one wavelength 
should be close to 760 nm while the other one can be within the range 830-900 
nm without substantially affecting the condition number. 

20 Using the spectra shown in Fig.1, model tissue absorption spectra were 
generated. Based on matrix inversion values of [oxy-Hb], [deoxy-Hb] and S were 
backcalculaled and the sensitivity to incorrect assumptions about the [water] and 
[lipid] tested. One approach is to take the measured na spectra and subtract 
water and lipid absorption corresponding to certain assumed concentrations. For 

25 the data shown in Fig. 4, a model tissue containing 15 jiM [HbTj, (true) 
saturation values of S =• 25%, 50% and 75% was used. Lipid concentration was 
40%, It was tested how a misjudgement of water concentration affects the 
recalculated S value. To test the error in a simple two-wavelengths-fit (760 and 
850 nm), the assumed lipid concentration was varied between 0 and 100%. 

30 When the assumed water concentration is right (lower three lines in Fig. 4), the 
deviation in saturation between true and calculated values is < ±2% (obviously 
with zero error for the right lipid concentration of 40%). A misjudgement about 
the water concentration by 20% (upper lines in Fig. 5) results in additional errors 
in S(calc)-S(true) of up to 2% for S=75%, 4% for S=50% and 8% for S=25%. 

35 These errors in S are a function of the underlying tissue absorption coefficients. 
The values here giv an indication about the order of magnitude. 
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Having a system with more than two fit-parameters, best wavelength 
combinations, for a three-components system of oxy-Hb, deoxy-Hb and lipid 
system, for a four-components system of oxy-Hb, deoxy-Hb f lipid and water, and 
the sensitivity of calculation of S to noise at the different wavelengths have been 
5 determined. 

in Figs. 5A and 5B the inverse of the condition number is plotted for a three 
wavelengths system based on the oxy-Hb, deoxy-Hb and Jipid specific 
absorption spectra as a function of X<i and X 2 . The third wavelength was fixed at 
73 = 830 nrn. Again, the plot is symmetric with respect to the diagonal. From Fig. 
10 5A it is apparent, that there are three "islands" of high 1/C values. Unfortunately, 
all of these island would include wavelengths outside an imposed hardware 
constrained wavelength range of 750 to 850 nm. Plotting the same data in a 
different scale (Fig. 5B) shows that there is just a single preferential combination 
within this hardware constrained wavelength range; 760 and 780 nm. 

15 Equivalent to Figs. 5A and 5B, the inverse of C for a 4-waveIengths system is 
plotted In Fig. 6A and 6B. Again, the difference between them is the scaling. 
Two wavelengths were fixed at 7J3 = 760 nm and M = 830 nm. Including the 
wavelengths outside the 750-850 nm range there appear four preferential 
combinations. Restricting the wavelength range to 750-850 nm there are just two 

20 advantageous region (marked by the white rectangle in Fig- 6B): 780 nm and 
850 nm, and 780 nm and 815 nm. 



From the analysis based on matrix condition numbers, the best wavelength 
combinations for 2, 3 arid 4 wavelengths measurements are the following: 
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Table 1 

Wavelength combinations for 2, 3 and 4 wavelengths measurements 





Wavelength 


Best wavelengths (nm) 




fit for 


see 




Range 


X1 


X2. 




X4 




Fig. 


2-a. 


650-950 nm 


700 


>860 






0*y-Hb, deoxy^Hb 


2.3 




750-850 nrn 


760 


850 










3-X 


650-950 nm 


700-760 


830 


925 




+ lipid 


5A.B 






830 


860-870 


525 










750-650 nm 


760 • 


780 


830 








A-\ 


650-950 nm 


760 


830 


860 


925 | 


+ lipid, + water 


6A,B 






700 


760 


830 


925 








750-650 nm 


760 


780 


830 " 


850 


Host combination 








760 


780 


815 


830 







5 Furthermore, it must be pointed out that including more wavelengths does not 
increase the condition number. Eg. for the four chromophores and all 
wavelengths in the range 750-850 nm, 1/C - 0.000314. This is lower than the 
value (1/C = 0.00036, compare with Fig. 6B) when only four wavelengths (760, 
780, 830 and 850 nm) are used. In a system without noise and no other 
10 chromophores than the four considered here, a 4-wavelengths system is the 
optimal. 

While certainly only a 4-wavelengths measurement allows [oxy-Hb], [deoxy-Hb], 
TJipid] and [water] to be determined, and a 2-wavelengths system (see Fig. 5) is 
not sufficient, the question is posed whether a a^wavelengths measurement 

15 might supply S values with a high enough precision. In this case the 
concentration of one chramophore (water or lipid) must be guessed and the 
conesponding absorption subtracted from the measured m values. This was 
tested with a model absorption spectrum and is shown in Fig. 7 for 760, 780 and 
850 nm. True water concentration was varied between 10 and 100% (the 

20 different lines), and the difference between calculated and true saturation values 
plotted as a function of assumed water concentration. For instance, for a true 
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water concentration of 50%, a misjudgment of the water concentration by 10% 
results in an error in S by about 4%. 

Up to now, only "perfect" data sets were considered with no noise. In real 
situations there are problems due to measurement noise that is random for the 
5 different wavelengths; unknown chromophores in the tissue, i.e. there Is a 
background absorption coefficient the spectrum of which we do not know; and 
possible systematic errors in the primary recovery. 

There are an ample number of parameters which can be considered and as 
examples the following two questions are considered. First, is the oxygen 
10 saturation more susceptible to noise at certain wavelengths? Second, what is 
the influence of an offset in the jia data? 

The influence of errors (noise) in on the calculated Hb concentrations and 
saturation values was estimated in a model tissue based on 20fiM [HbT], S= 
50% and a lipid and water concentration of 30% and 40% respectively. Matrix 

15 Inversion was performed on the ^ - spectrum of this model tissue for 
wavelengths 760 f 780, 830 and 850 nm. The sensitivity to noise (i.e. variations 
in at the different wavelengths was estimated by changing the absorption 
coefficient at a single wavelength by +0.0001 mm" 1 . In Fig. 8, the change in 
calculated [oxy-Hb], [deoxy-Hb] and oxygen saturation value due to this "noise" 

20 is plotted. This figure shows that the change in oxygen saturation value is about 
-2% for changes at 760 nm, <0,5% at 780nm, while it translates to a variation of 
+6% at 830 nm. 

There are two conditions that produce an offset in the measured jia-spectra with 
respect to the true values. First, the algorithm for ^calculation based on TPSF- 

25 based optical imaging might lead to a systematic offset e.g. due to residual 
crosstalk between absorption and scattering parameters. Second, the tissue 
absorption might have a background of unknown origin (chromophore). Under 
both conditions the fitting of j%-data with the four chromophores is hampered. 
The effect of such an offset for different wavelength combinations is estimated 

30 with a model spectrum of 20 \iM [HbT], S=75% and water and lipid concentration 
of 40%. An offset of 0.0005 mm -1 was added to the na-values independent of 
wavelength. The effect on the calculated oxygen saturation values is known in 
Fig. 9 for combination of 2, 3 and 4 wavelengths as well as continuous spectra 
between 750-850, 720-850 and 720-900 nm. It Is apparent that the lowest error 

35 in S is achieved by the 4-wavelengths combination. Including more wavelengths 
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increases the error. In Fig. 10 the same calculation was done, however, for a 
true oxygen saturation value of 50%. Here the lowest error is achieved by the 
720-850 nm wavelength range, while using less wavelengths or increasing the 
fitting range to 900 nm results in larger errors. 

5 Based on the assumption that the dominant tissue chromophores are oxy-Hb, 
deoxy-Hb, water and lipid and analysis of the matrix condition number, 
measurements at the wavelengths 760, 780, 830 and 850 nm supply an optimal 
data set when the wavelength range is limited to 750-850 nm under Ideal 
conditions. As shown in Figs. 5 and 6, inclusion of shorter and longer 

10 wavelengths promise a better matrix inversion. Under real conditions there is no 
clear-cut answer about the improvement when more wavelengths are included 
(see Figs. 9 and 10). It might be advantageous to reduce measurement noise at 
4 wavelengths due to longer scan times rather than to include more 
wavelengths. As demonstrated in Fig. 8, to achieve an optimal accuracy the 

15 noise level at different wavelengths has to be adjusted which might require 
different measurement times at certain wavelengths. 

Strictly speaking the work presented here was achieved by optimizing a 2- 
wavelength system and then optimizing a 4-wavelength system where 2 of the 
wavelengths were fixed at the optimized 2-wavelength system values. Whilst this 
20 is easier to display graphically, preferably all 4 wavelengths would be permitted 
to vary in a global optimization process. Fortunately, for the specific example 
presented here when all 4 wavelengths are permitted to vary the same optimal 
solution is found. However, this may not be true for all situations and a global 
optimization is preferred. 

25 It is also understood that it will be obvious to those skilled in the art that the 
same approach for choosing optimal wavelengths can be applied to optical 
absorption spectroscopy in general. For example, in other embodiments of the 
present invention the method of the present invention is also used for choosing 
the optimal wavelengths for analyzing the components of paints, pharmaceutical 

30 products, food, grain or any other turbid media. 

It is also understood that the proposed method applies both to the analysis of 
* absolute chromophore concentrations as to their changes or relative 
concentrations. 

While the invention has been described in connection with specific embodiments 
35 thereof, it will be understood mat it is capable of further modifications and this 
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application is intended to cover any variations, uses, or adaptations of the 
invention following, In general, the principles of the invention and Including such 
departures from the present disclosures as come within known or customary 
practice within the art to which the Invention pertains and as may be applied to 
5 the essential features herein before set forth, and* as follows in the scope of the 
appended claims. 
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CLAIWS 

1. A method of optical imaging of turbid media using a plurality of 
discrete wavelengths in an optical imaging system, the method comprising the 
steps of: 

selecting a set of chromophores for characterizing a property of the 
turbid media; 

defining parameters of the system including at least a number of said 
discrete wavelengths, a value of each of said wavelengths, source power and 
detector aperture for each of said wavelengths, a choice of image algorithm and 
source/detector geometries, a choice of source and detector and noise 
characteristics; 

fixing a value of all of said parameters except a plurality of said 
parameters values to be optimized; 

determining an optimal value for each of said parameter values to be 
optimized as a function of a performance of the system in measuring a 
concentration of said chromophores in said turbid media for characterizing said 
property as a whole; and 

using said optimal value for each of said parameter values in imaging 
said turbid media. 

2. The method of claim 1, wherein said imaging is medical imaging, said 
highly turbid medium being body tissue and said property is physiological. 

3. The method of claim 1 or 2, wherein said parameter values to be 
optimized comprise a value of each of said wavelengths. 

4. The method of claim 3, wherein said parameter values to be optimized 
further comprise said number of said discrete wavelengths. 

5. The method of claim 4, wherein said step of determining comprises 
fixing said number of discrete wavelengths at each of a plurality of numbers, and 
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determining an optimized performance of the system in measuring a 
concentration of said chromophobes in said turbid media at each of said plurality 
of wavelengths, and selecting one of said plurality of numbers having a best 
optimized performance. 

6. The method of claim 3, wherein said step of determining an optimal 
value for each of said parameters comprising minimizing a condition number of 
a matrix of specific absorption coefficients of said chromophores as a function of 
wavelength- 

7. The method of claim 4, wherein said step of determining an optimal 
value for each of said parameters comprises minimizing a condition number of a 
matrix of specific absorption coefficients of said chromophores as a function of 
wavelength. 

8. The method of claim 5, wherein said step of determining an optimal 
value for each of said parameters comprises minimizing a condition number of a 
matrix of specific absorption coefficients of said chromophores as a function of 
wavelength. 

9. The method of claim 1, wherein said step of determining comprises 
empirically detemnining said performance of the system for a range of said 
values for each of said parameter values to be optimized.- 

10. The method of claim 2, wherein said plurality of chromophores 
comprise at least oxy-hemoglobin and deaxy-hemoglobin. 

11. The method of claim 10, wherein said chromophores are water, lipids, 
oxy-hemoglobin and deoxy-hemoglobin. 



12- The method of claim 1 0, wherein said body tissue is breast tissue. 
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13. The method of claim 10, wherein said number of wavelengths 
selected is from 2 to 4. 

14. The method of claim 1 3, wherein said number Is 4. 

15. The method of daim 11, wherein values of said wavelengths are 760 
nm, 780 nm, 830 nm and 850 nm. 

16. The method of daim 1 t wherein the step of determining an optimal 
value of said parameters to be optimized comprises: 

deriving an inherent wavelength-dependent sensitivity to noise in 
calculating said chromophora concentrations, and 

determining an optimal correlation of said sensitivity and at least one 
other of said parameters. 

17. The method of claim 16, wherein one of said parameters to be 
optimized is a distribution of an acquisition time at each of said wavelengths. 

18. The method of claim 1, wherein one of said parameters to be 
optimized is a distribution of an acquisition time at each of said wavelengths. 

19. The method of claim 18, further comprising a step of determining a 
minimum value for said acquisition time at which said performance of said 
system attains a minimum threshold value. 

20. The method of claim 1, wherein one of said parameters to be 
optimized is at least one of said source power and said detector aperture for 
each of said wavelengths. 
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21. The method of claim 20, further comprising a step of determining a 
minimum value for an acquisition time at which said performance of said system 
attains a minimum threshold value. 
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Fiber-optic probe for noninvasive real-time det rmination 
of tissue optical properties at multiple wavelengths 



Jan S. Dam, Carsten B. Pedersen, Torben Dalgaard, Paul Erik Fabricius, 
Prakasa Aruna, and Stefan Andersson-Engels 



We present a compact, fast, and versatile floor-optic probe system Tor real-time determination of tissue 
optical properties from spatially resolved continuou»wave diffuse reflectance measurements. The ays- 
lea collects one act of reSectance data from six source-detector distances at four arbitrary wavelengths 
with a maximum overall sampling rate or 100 Hz. Multivariate calibration techniques based on two- 
dimensional polynomial fitting: are employed bo extract and d lap lay the absorption and reduced scattering 
coefficients in real-time mode The foui wavelengths of the current configuration are 660. 785, 805. and 
974 nm. respectively. Cross-validation teals on a 6 x 7 calibration matrix of IoJfcralipid-dye phantoms 
showed that che mean prediction error at. e.£ M 786 nm was 2.8% for the absorption coefficient and U* 
for the reduced scattering coefficient The errors are relative to the range of the optica] properties of the 
phantoms at 78S nm, which were 0-0-3/cm for the absorption coefficient and 6-16/cm for the reduced 
scattering: coefficient Finally, we also present And discuss results from preliminary skin Ussue mea- 
surements. © 2001 Optical Society of America 

OCIS codes: 170.1470, 170.1580. 170.3890, 17OJ6510. 



1. Introduction 

The optica] properties of human tissue, 2 i.e., the ab- 
sorption coefficient tbe scattering coefficient 
and the alnisotropy factor g t may provide important 
information on the composition and the physiological 
dynamics! of the tissue. Whereas u. c may provide 
information on tissue chr omophores, 3 - * and g may 
be used to characterize the form, size, and concentra- 
tion of various scattering components in the cis- 
sue. 5 ~ a | 

Owing |to the obvious advantages of noninvasive 
and minimally invasive measurements, determina- 
tion of tissue optica] properties based on diffuse re- 
flectance measurements has a significant potential in 
biomedical diagnostics and monitoring Diffuse re- 
flectance measurements may be roughly divided into 
time-resolved 9 - 10 (TR), frequency-domain 11 ^ (FD), 
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and spatially resolved continuous-wave 1 *"* 1 (cw) 
methods. Traditionally, TR- and FD-based methods 
have been considered to be more accurate for absolute 
determination of optical properties than cw based 
methods. However, TR and FD methods also re- 
quire more bulky and expensive equipment, and a 
larger sample volume, which may restrict some bio- 
medical applications, e.g, p implementation in porta- 
ble monitoring equipment or in endoscopes and 
catheters. Recent research 19 has shown that cw- 
based methods may yield absolute determination of 
the optical properties of tissue with accuracies siini. 
lar to the TR and FD methods, which makes cw -based 
methods a better choice for many practical applica- 
tions. 

Continuous-wave diffuse reflectance methods may 
be further divided into methods based on probes in 
contact with the tissue and noncontact methods, i.e. t 
image rcflectometry. The latter method is advanta- 
geous in clinical applications because of the noncon- 
tact and thus sterile properties. The advantage of 
contact probes is that they can be made small and 
portable. Therefore they are well suited for (a) op- 
tical biopsies of body cavities or organs, because a 
probe may be implemented in existing endoscopic 
equipment, or (b) for long-time monitoring of tissue 
optical properties, because tba prob may be fixed to 
the skin and still allow the patient to move around. 
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Fig. 1. Diagram of the fiber-optic system for Rir) measure on ante applied in this paper. (I) Probe head with source and detector optical 
fibers mounted in a rotational symmetric configuration. <2) Handheld box iviih silicon photodiodes and amplifier electronics. (33 
Stationary box con lain in? a distal signal processing (DSP) board and Che light sources In the form of diode lasers. (4) External 
temperature controller to maintain a constant temperature of the diode lasers. (5 1 Laptop PC to analyze, display, and store the acquired 
Kir) data. 



However, there arc a number of drawbacks related to 
skin tissue monitoring using" a contact probe system. 
First, the static mechanical pressure of the probe 
may influence the tissue optical properties, e.g., by 
displacement of blood in the sample volume. Sec- 
ond, physical activities of the patient may lead to 
motion artifacts, Le., movements of the probe relative 
to the skill tissue. Finally, obstruction of the tissue 
and sweating may cause long-term changes in the 
tissue optical properties and the probe^skin optical 
interface. All such contributions have to be isolated 
to extract the relevant information from diffuse re- 
flectance measurements. 

Accurate closed-form mathematical analytical ex- 
pressions for spatially resolved diffuse cw reflectance 
2?(r) as a function of the optical properties axe 
strongly limited by requirements to the range of the 
optical properties and the specific geometry of the 
setup. 22 Therefore methods based on more-accurate 
numerical light-propagation models, e.g., Monte 
Carlo simulations, 23 are often used in conjunction 
with various forms of multivariate analysis 2 ** 24 to 
solve the inverse problem of extracting the optical 
properties from R(r) measurements. However, as a 
consequence of the similarity principle, 2 * 2 * normally 
only and the reduced scattering coefficient - 
Cl - are extracted from R(r) measurements. 



In this paper we present a fiber-optic probe system 
for in vivo real-time determination of tissue optica] 
properties based on R(r) measurements at four arbi- 
trary wavelengths. The probe system collects the 
diffuse reflectance at six distances, extracts the cor- 
responding and u/, and displays the results im- 
mediately on a laptop personal computer' (PC). The 
prediction of and \x. s ' is based on multiple polyno- 
mial regression 27 and calibration on a set of 
Intralipid-dye phantoms with well-known optical 
properties within a range typical for biologic tissue 28 
in the visible and the near-infrared regions. In the 
fallowing sections we first give a description of the 
probe specifications and the principles of the calibra- 
tion and prediction algorithms. Next, we present 
and discuss the results obtained from phantom mea- 
surements and simulated numerical tests with re- 
spect to prediction errors, robustness, and the like. 
Finally, we present the results of preliminary predic- 
tion tests obtained from a series of clinical measure- 
ments on normal skin tissue. 

2. Materials end Methods 

A Rber Probe System 

Figure 1 shows the basics of the fiber probe system we 
used for the measurements described below. The 
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system consists of a probe head with a source 
fiber in the center surrounded by five equally spaced 
concentric rings of 250-iim detector fibers. We chose 
this ring geometry instead of, e.g., a simpler linear 
geometric configuration, in part to be able to collect 
more light at farther distances and in part to mini- 
mize any problems arising from tissue inhomogene- 
ifcies during clinical measurements. The fibers of 
each single ring detector are bundled and terminated 
on separate silicon photodiodes. In addition, three 
photodiodes and e temperature sensor are mounted 
directly near the perimeter of the probe head- Thus 
J?(r) can be collected at six distances, in part at the 
five fiber rings (r « 0.6, 1.2, 1.8, 2.4, and 3.0 mm) and 
in part at the directly .mounted photodiodes (r =■ 7.8 
mm). These distances were chosen on the basis of 
previous studies, 19 The gain of each reflectance de- 
tector was calibrated in an integrating sphere setup 
to obtain equal outputs at constant input light inten- 
sities- The source fiber is coupled into four separate 
fibers, each connected to four replaceable low-power 
diode lasers. The diode lasers are mounted on a 
heat sank with a constant temperature maintained by 
an external controller. Furthermore, a separate ref- 
erence detector monitors the output of the source 
fiber at the probe head. The diode lasers may be 
selectjed arbitrarily to suit different applications. In 
this paper we used diode lasers with the wavelengths 
660, 785, 805, and 974 urn, which are well suited for 
applications involving hemodynamic monitoring. 
The data acquisition and storage is controlled by a 
laptop PC connected to a digital signal processing 
(DSP) board. In each R(r) measurement the detec- 
tor hardware collects data simultaneously in eight 
parallel channels from the probe head, i.e., from the 
six detector rings, from the reference detector at the 
source fiber, and from the temperature sensor. One 
cycle of four successive measurements (i.e., one at 
each wavelength) including dark measurements may 
be performed in —10 ms; thus the maximum sam- 
pling rate of the system is —100 Hz. To minimize 
any interference from background light or drift of the 
light source, the dark measurements are subtracted 
from the measured reflectance data after which they 
are normalized relative to the source reference. The 
DSP board accomplishes this prior to when the data 
are analyzed, displayed, and stored by the PC. 



B. Calibration and Prediction Algorithms 
In theory, and may be determined with R(r) 
data from only two of the six detector distances of the 
fiber probe. Building on our previous research, 1927 
we thus applied multiple polynomial regression 
(MFR) to create a calibration model and subsequently 
extract ^ and \*J from R(r) measurements atr, = 0.6 
mm and at r z = 7.8 mm. for the sake of clarity we 
give a summary of the MPS method here. We first 
measure J?(r) at rj and at r 2 for a set of calibration 
samples with well-defined optical properties and de- 



note them R ltCo] and R 2xa) - Then we find a double- 
polynomial fit to J? liM ) and R^c&v 

x (rf fl + dm,/ 2 + . , ■ <W), O) 

Where a, 6, c, and d are fitting coefficients deter- 
mined by leasUsquares regression. J?i* L **d R^n 
thus constitute the calibration model. The next step 
is to solve the inverse problem of determining p tf and 
jjl,' from R(r) measurements on a set of prediction 
samples, i.e., R^ mas and J? 2 ,mo*r First we define 



(2) 



Then we use a Newton-Raphson algorithm to per- 
form converging iterative calculations of and p./; 



[F(n 0>> m/)V 



3F 3F 
BG dC 



te) 



{::::) -fcWd J 



(3) 



where k is an index and h Q and are correction 
terms of p. a and The iterations were stopped 
when both h 0 and h s < 1 X 10"* (cm -1 ). With such 
strict criteria, previous research 27 showed that the 
contribution of the Newton-Raphson method to the 
total prediction error was negligible. 

Since in this study we want to extract two optical 
properties only, i.e., ^ and ji/, the MPR method 
implies exactly two input variables as well, i.e., R(rJ 
and R{r 2 ). However, owing to tissue inhomogeneity 
and/or noisy measurement conditions, it might be 
advantageous to include more source-detector dis- 
tances daring measurements and then subsequently 
apply some sort of dimension-reduction method be- 
fore the data are fed into the MPR method. To test 
such a procedure, we also applied principal compo- 
nent analysis (PCA) on the data from all six source- 
detector distances of the probe system and then used 
the resulting two main principal components and 
jP 2 as input to the MPR method instead of iffo) and 
R(rJ. m 

All the prediction algorithms we applied in this 
paper were implemented in Matlab and run on a 
16&-MHz PC. With this configuration, the predic- 
tion of a single set of m, and could be performed in 
~60 ms. 
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Tabfo 1. Opltefll Property Ranges of S*7 Matrix of intmOpUWnk 
PnonJcm* Determined from Integral* g Spher* Moasuroroento 





Optical Properly Range (em" 1 ) 


X(nm) 






6S0 


0-O-36 


7.3-19.5 


785 


0-0.32 


8.1-18.3 


BOG 


0-0.31 


6.0-16.0 


974 


0.45-0.68 


4.S-12.7 



C. Phantoms 

Because of the unknown numerical apertures of the 
fiber probe light source and detectors, we chose to 
calibrate the system directly on a set of phantoms 
instead of using a mathematical light-propagation 
modeL The phantoms consisted of well-defined 
aqueous solutions of Intralipid and black ink in cy- 
lindrical glass containers with a diameter of 10 cm 
and a height o!4 cm, We determined the scattering 
and absorption spectra of the Intralipid and the black 
ink from integrating sphere 27 and traditional trans- 
mission spectroscopy measurements. On the basis 
of these spectra, we mixed a 6 x 7 matrix of phantoms 
with p*, and ranges typical for skin tissue (see 
Table 1). The applied ranges of Intralipid concen- 
trations were 0.6, (18, . . . , 1.6%, and the range of the 
ink concentrations were 0.0. 0.2, . , . , L29k It 
should be noted that the absorption of pure ink is 
much higher than Oiat of typical biologic substances; 
thus the ink concentrations in the following refers to 
a premixed basic ink-water solution with a biologi- 
cally relevant absorption level. To perform the pre- 
diction experiments, we calibrated the probe system 
directly to the concentrations of the Intralipid and 
the ink in the phantoms, assuming that both the 
absorption of pure Intralipid 2 *' 80 and the scattering of 
the ink were negligible. At 660, 785, and 805 nm we 
assumed that the background absorption was negli- 
gible as well. However, water exhibits substantial 
absorption 31 at 974 nm (m, fWMaT = 0.45 cm" 1 }. This 
was incorporated in the concentration-to-absorption 
tabulation by assuming a constant water absorption 
contribution at this particular wavelength (i.e.. = 
0%=^^ = 0.46 cm" 1 ). 

3. Results 

Once the 6 x 7 matrix of phantoms was mixed, we 
determined the actual optical properties of a subset of 
the phantom matrix, using integrating sphere mea- 
surements. Figure 2(a) shows the mean absorption 
spectra of the six phantoms with 1% ink concentra- 
tions but varying Intralipid concentrations, and Fig. 
2(b) shows the mean scattering spectra of the seven 
phantoms with 1% Intralipid concentrations but 
varying ink concentrations. We used the results 
from Fig. 2(a) to relate the concentrations of In- 
tralipid and ink to the optical properties of the phan- 
toms. The resulting optical property ranges of the 
6X7 matrix phantoms are listed in Table 1. 
Figures 3(a) and 3(b) show the measured intensity 
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Kg. 2. Optical property spectra of the applied Infcralipid-iok 
phantoms, (a) Mean of as a fundi on or the wavelength X lor 
the six phantoms with l<£ ink concentration but varying Intralipid 
concentrations, (b) Mean of m* as a Junction of the wavelength for 
the «even phantoms with 1* Intralipid concentration bat varying 
ink concentrations. The error ban of the spectra in panels (a) and 
(b) «bow the standard deviations of m, and u«\ respectively. The 
circles indicate the values at the four wavelength* of the probe. 
The fevo spectra were measured and calculated with an integrating 
sphere setup in conjunction with the MPfl method. 



R{r) at 785 nm as a function of ji ? and ft,' at r = 0.6 
mm and r = 7.8 mm, respectively. Corresponding 
plots at the three remaining wavelengths of the 
probe, Le., 660, 805, and 974 nm, showed similar 
characteristics; ie., they were also smooth and mono- 
tonic, which suggests that they may be fitted well by 
Eq. (1) with relatively low-order polynomials. Pre- 
liminary experiments showed that using Eq. (1) with 
m - 3 provided the best overall calibration model, 
concerning accuracy and robustness. The prediction 
accuracy of the MPR method was tested with leave- 
one-out cross validation. This means that we suc- 
cessively performed predictions with the data from 
one phantom for prediction and the data from the 
remaining- 41 phantoms of the 6 x 7 matrix for cali- 
bration. To ensure that all R (r) variations were cov- 
ered by the calibration models, we carried out 
predictions only on the 4X5 interior subset of the 6 X 
7 matrix. Table 2 presents the results from these 
cross-validation prediction tests. 
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(a) Phantom mBasuremenfc al r= 0.5 mm W Phantom measurements at r= 7 J mm 




(c> stmulaied valuee at r* 0.8 mm < d > Simulated values at r = 7, a mm 




Tig. 3 Surface plots olRir) at r x = (L6 mm and r z = 7.8 mm as a Jbnction of m, and >0 Panels la) and (b) show the R{r) plots of tite 
6 * 7 Inindipid phanioms at 78S am, and panels <c) and (d) show tha corraspondine *(D plots based on Monte Carlo simulabona. Note 
that the arbitrary inleasJiy units of measured R(r) plots have been scaled in order to compare them with the simulated plots. 



To check our experimental results, and to examine 
the robustness, noise sensitivity, and so on, of the 
applied calibration model and prediction algorithms, 
we also performed a series of numerical tests based 
on Monte Carlo simulations with a geometric config- 
uration analogous to the fiber probe setup described 
in Subsection 3. A. However, unlike the measure- 
ments, the simulations were baaed on a collimated 
incident beam, a semi-infinite medium, and 1 /nr ac- 
ceptance angle far collection of the diffuse reflectance. 
Hie tests were carried out as follows: (I) First, we 
generated a basic 6x7 calibration matrix with op- 
tical properties corresponding to the 6 x 7 phantom 

Table 2. Leave-One- Out Cmsa-VoGdaUon Prediction Tests Baaed on 



Phantom Mtasura menu" 







Predicti an E/ror* ('* \ 




Pa 








X(nm) 


Mean 


Max 


Meon 


Max 


600 


3.3 (2.5) 


18 (6.7) 


1.7fl"i 


.1.714.01 


785 


2A 


8.3 
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2A 


BOS 


2-6 


9.0 


1.5 


AS 


974 


3.7 


9.7 


1.6 


n.K 



"The errore are relative to the ranges of m and a/ of the phnn- 
turns. The bracketed values at X ° 660 nro beva been corrvi-hiJ for 
aa outlier detected by visual inspection of too measured Hir\ HnUi. 



matrix at \ = 785 nm and then performed prediction 
tests similar to the phantom cross-validation tests 
described above. Figure 4 shows a selected set of 
R(r) profiles from these simulations [see also Figs, 
3(c) and 3(d)]. (II) Next, to teat the algorithms with 
independent calibration and prediction data, we gen- 
erated a series of 20 R{r) data sets with random ^ 
and jij/ distributions and then performed predictions 
on this randomized set with a calibration model 
based on the full 6 X 7 basic matrix." (TH) The basic 
6X7 matrix was generated with 1 x 10* photons. 
To examine the impact of random measurement noise 
on the prediction performance, we generated a second 
6X7 matrix with less noise by using 1 x 10 7 photons 
and performed identical cross-validation tests on this 
set (TV) We also wanted to investigate the effect of 
using a calibration model with a higher u*, and ^ 
resolution; thus we also carried out cross-validation 
tests on a II x 13 calibration matrix with the same 
u 4 and range and number of photons as the basic 
6X7 matrix. (V) Finally, we tested the effect of 
using PCA for dimension reduction as described in 
Subsection 3.B. The resulting mean and maximum 
prediction errors from the numerical tests Q-V) de- 
scribed above are summarized in Fig. 5. 

To validate the optical property range of the ap- 
plied calibration model, we performed a series of 
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Fig. 4. Monte Carto simulated H{r) data for various wrnhinaticmB cf ^ and ^' within ranges typical for skin tissue at 785 ran. In Ul 
and (b), is kept constant at 6 cm" 1 and 16 cm~\ raspectively while is varied within th« range 0-0.3 an" 1 . In to a&d (d) is VcpL 
«Mtant at 0 em"* and 0.3 can" 1 , respectively, uhercas is varied within the range 6-16 cm" 1 . 



probe measurements on the inner forearm of Ave 
healthy individuals. Figure 6 shows how the mean 
R(r) of all five individuals (rircies) at each probe dis- 
tance is situated relative to the J?(r) ranges of the 
calibration model (vertical bars). Figure 6 also 
shows the predicted optical properties at all four 
wavelengths calculated from the mean J?(r) data with 
Eq. (l)-(3). It can be seen that the measured skin 
tissue at X = 660 nra and r = 7.8 mm is outside the 
calibration model range, and therefore it is not pos- 
sible to predict \i a and p./ correctly in this case. 

4. Discussion 

A Geometry Considerations 

In Fig. 4(a) end 4(b) p*' is kept constant at values of 
6 and 16 cm" 1 , respectively. In these two cases, it 
appears that changes in ^ have a negligible effect on 
R(r) at distances close to the source. In Figs. 4(c) 
and 4(d), p* is kept constant at values of 0 and 0,3 
cm" 1 , respectively. Here it is notable that there is 
very litUe variation in i*(r) at r ** 0.35 cm in Kg. 4(c) 
and at r ~ 0.3 in Fig. 4(d). Using a figure of speech, 
we say that there is a pivot point in the R(r) graphs at 
r « 0,3 cm, when jj^ is kept constant. The simula- 
tions in Figs. 4(a) and 4(b) indicates that p./ may be 
determined with good accuracy from dose distance 
measurements only. To determine p* as well, Figs. 
4(c) and 4(d) suggest that R(r) measurements close to 
the pivot point should be included, since there is little 



variation in 2?(r) as a function of m' at this point. 
Although, other authors 32 - 33 support this argument, 
we based our experiments in this paper on dose 
range distances in conjunction with distances well 
beyond the pivot point, We did this because our 
previous studies 18 showed that this geometrical con- 
figuration provided a better accuracy than a config- 
uration with close range distances in conjunction 
with distances near the pivot point. 

B. Phantom Measurements 

The prediction tests using Intralipid-ink phantoms 
show a good accuracy (see Table 2). The mean pre- 
diction errors at all four wavelengths is roughly 3% 
for m,, and 1.6% for m\ whereas the maximum pre- 
diction errors for ^ a and are approximately 11% 
and 4%, respectively. The prediction algorithm con- 
verged in all cases; however, the maximum error of \l c 
at X = 660 ran is remarkably high (18%). A visual 
inspection of the raw J?(r) data at 660 nm revealed an 
outlier at r - 7.8 mm in one of the phantom mea- 
surements. This outlier could be due to an air bub- 
ble or a piece of dirt at one of the three 7.8-xom 
detectors in this particular measurement When we 
removed the outlier from the calibration and predic- 
tion analysis at X = 660 nm, we obtained the im- 
proved results shown in brackets in Table 2. In 
general, the prediction errors of p^ ore approximately 
twice as high as the errors of \i t \ This may be at- 
tributed in part to the fact that \i a is mainly deter- 
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rained on the basis of the j?(r) data at r =» 7.8 mm, 
whereas m' is almost solely determined from R{r) 
data atr = 0-6 mm, where the signal level is —1000 
times the level at r = 7.8 mm- The ^ predictions 
are therefore more sensitive to any background noise 
interference during the measurements. 

With reference to the comments on the pivot point 
in Subsection 5.A we also performed a series of pre- 
diction tests, in which we replaced the R(r) probe data 
a t r z = 7.8 mm with probe data at the pivot point r piwot 
= 8.0 mm. This had no significant effect on the p./ 
prediction errors, but both the mean and the maxi- 
mum jj^ prediction errors increased —30% compared 
with results reported in Table 2. 

C. Numerical Tests 

The general shape of the measured and the simulated 
R(p. a , plots in Fig. 3 are very similar at corre- 
sponding detector distances. However, the simu- 
lated plots do show slightly higher #(ji a> m*) levels 
for low and \l s ' values at r = 7.8 mm. This dif- 
ference is severe enough to prevent the direct use of 
a calibration model based on Monte Carlo simula- 
tions for prediction analyses of R[r) data collected 
with the fiber probe. The difference between the 
simulated and the measured R{r) data in Figs. 3(b) 
and 3(d) may be attributed to the different numerical 
apertures of the sources and the detectors in the two 
cases; e.g., the Monte Carlo simulations employ a 



collimated source beam, whereas the probe source 
beam is divergent. 

The r> 0 prediction errors of the simulated data 
shown in Figs. 5(a) and 5(c) are comparable with the 
errors of the measured data Table 2, whereas the >a/ 
prediction errors of the simulated data [Figs. 5(b) and 
5(d)] are lower than the errors of the measured data. 
In both cases the \i a prediction errors are signifi- 
cantly higher than the m/ errors, probably owing to 
the noise sensitivity at r = 7.8 mm as discussed 
above. Figure 5 also shows a drop in the prediction 
errors when a randomized prediction data set is used. 
This may be because the quality of the calibration 
model fit during leave-one-out cross validation is re- 
duced in the vicinity of the ji a and p.,' to be predicted. 

From Fig. 5 it also appears that a reduction in 
measurement noise (Le., more photons) or use of a 
calibration model with a higher jl c and p*' resolution 
in both cases lead to a significant drop in maximum 
prediction errors, whereas the effect on the mean 
prediction errors is moderate. Finally, it can be seen 
from Kg. 5 that the use of PCA before calibration and 
prediction has a positive effect on the ^ prediction 
errorB, whereas the \jl 0 ' prediction errors are slightly 
increased. In summary, the prediction accuracies 
obtained from calibration and prediction on phan- 
toms are comparable with the accuracies obtained 
from calibration and prediction on Mimte Carlo sim- 
ulated data. The numerical tests suggest that the 
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prediction accuracy may be improved slightly by use 
of more samples during creation of the calibration 
model. However, the improvements may hardly jus- 
tify the additional efforts of creating more calibration 
phantoms. 

D. CRnical Measurements 

From Fig. 6 it appears that the range of calibration 
model covers the span of the measured tissue R(r) 
data in most cases, except at X = 660 nm. At this 
wavelength the mean of the measured R(r) at r 2 «= 7-8 
mm is above the range of the calibration model. 
This suggests that either n^arji,' (or both) of the skin 
tissue is lower than the minimum values of the cali- 
bration model. In the ca6e of p a this is not possible, 
since the model range at X = 660 nm is 0-0.36 cm -1 . 
The minimum value of |i/ of the model at this wave- 
length is 7.S cm" 1 , and it is not likely that the mea- 
sured p., should be beneath this limit either. At the 
three remaining wavelengths the measured J?(r) at r x 
= 0.6 mm and r 2 — 7.8 ™™ are within the limits of the 
calibration mode, and the extracted u ft oru,' shows a 
reasonable correlation with other studies. 20 ' 34 How- 
ever, at 974 nm, it appears that the measured R(r) 
data, at some of the intermediate distances (i.e., at 
r = 2.4 and 3.0 mm), are not within the ranges of the 
calibration model; Le,, the shape of the tissue R(r) 
profiles are different from the shape of phantom pro- 
files. There may be several reasons for this. First, 



in contrast to the phantoms, the volume of the skin 
tissue sampled by the probe is not homogeneous but 
basically consists of distinct layers, i.e., the stratum 
corneum, the vivid epidermis, the epidermis, and so 
on, Each of these layers has specific optical proper- 
ties, and therefore the predicted values of ^ or jjl/ in 
Fig, 6 are a compound of the optical properties of each 
separate layer. Second, the ratio of the refractive 
indices between the probe and the sample (n^bc/ 
*aampie) is different in the two cases; Le. t n damM ■* 1.40 
(Ret 35) and n phMlMO « 1.33. Finally, the probe- 
tissue interface may not have been optimal, which in 
turn may lead to light piping effects between the 
probe and the stratum corneum. Such effects may 
cause an increase of the farther R{r) values relative to 
the values close to the light source, which also might 
explain the out-of-range problems encountered at \ = 
660 nm. To summarize, it is essential to ensure a 
proper optical contact and apply a uniform and re- 
producible mechanical pressure in order to obtain 
valid predictions from contact probe measurements. 

E Acquisition end Analysis Speed 
As stated in Subsection 3, A, the maximum acquisi- 
tion rate of the fiber probe system for all four wave- 
lengths is —100 Hz, which is sufficient for analyzing 
the optica] property dynamics of most physiological 
systems. The current software implementation of 
the prediction algorithms leads to an analysis time 



1 162 APPUED OPTICS / Vol. 40, No. 7 f\ Marcn 2001 



01/16/2004 15:46 FAX 613 230 6706 



OGILVY RENAULT 



81112 



for one set of ^ or ^/ of -50 ms, i.e„ an analysis rate 
of 5 Hz when all four wavelengths are included To 
match the analysis rate to the acquisition rate of the 
probe system, the speed of the prediction algorithms 
may be increased substantially by application of a 
more contemporary PC and/or by compiling the al- 
gorithms in, e-g., the C programming language. 

a Conclusions 

We have demonstrated a versatile, fast, and accurate 
probe system for real-time noninvasive determina- 
tion of tissue optical properties from continuous wave 
(cw) spatially resolved diffuse reflectance measure- 
ments. The current calibration of the system was 
intended to match a typical range of absorption and 
reduced scattering coefficients of skin tissue, and pre- 
liminary clinical trials on a set of healthy individuals 
showed a good consistency except at X = 660 nro. At 
this wavelength the measured tissue data were out- 
side the range of the calibration model in some In- 
stances. This problem may be solved by extension of 
the calibration model range, but the problem may 
also originate in the fact that the system was cali- 
brated on homogeneous one-layer phantoms, 
whereas skin tissue is an inhomogeneous multilayer 
structure- Furthermore, the clinical trials clearly 
demonstrated that it is imperative to ensure proper 
and reproducible refractive-index matching and me- 
chanical contact properties to obtain valid results. 
All things considered, the system and the method we 
have presented here provide a sound basis for future 
development of compact and dedicated systems for 
noninvasive or minimally invasive medical diagnos- 
tics and monitoring. Still, further research is re- 
quired for exploring the applicability for specific 
biomedical implementations. 
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Stains have been used in opiical microscopy to visualize the distribution and intensity of substances 
to which they are attached, Quantitative measures of optical density in the microscopic images can 
in principle be used to determine the amount of the stain. When multiple dyes are used to simul- 
taneously visualize several substances lo which they are specifically attached, quantification of each 
stain cannot be made using any single wavelength because attenuation from the several stain 
components contributes to the total optical density. Although various dyes used as optical stains are 
perceived as specific colors, chey. in fad. have complex attenuation spectra. In this paper, we 
present a lechnique for multiple wavelength Image acquisition and spectral decomposition based 
upon the Lambert-Beer absorption law. This algorithm is implemented based on the different 
spectra) properties of the various stain components. By using images captured at N wavelengths, N 
components with different colors can be separated This algorithm is applied to microscopy images 
of doubly and triply labeled prostate tissue sections. Possible applications are discussed. © 1996 
American Association of Physicists in Medicine. 



I. INTRODUCTION 

In pathology, human tissue sections are commonly stained 
biochemically or immunccyiocherajcalfy in order to use 
color to enhance the contrast of specific structures in the 
microscopic images. The distribution of stain provides an 
indicator of the distribution of the substance or structure to 
which the stain is specifically attached. The optical density 
of the stain can then be used to determine the amount of 
stain. This measurement can be simply obtained through a 
monochrome image. 

In nearly every situation, pathologists desire to simulta- 
neously evaluate a number of substances or structures in a 
single piece of tissue. There are generally two purposes for 
the evaluation of multiply stained slides: (1) To evaluate a 
specific substance in the structure that can only be visually 
recognized through the second dye. For instance, in order to 
study the distribution of immunocyiochernicaJ PSAP staining 
in prostate tumor cells, a second nuclear stain is used for 
recognition of malignant cells. (2) To obtain the relationship 
of two or more substances or to simultaneously examine the 
photometric and morphomeuic parameters of one or more 
suruciures. For example, in the study or prostate cancer, the 
PSAP staining end nuclei shape may provide some diagnos- 
tic information. 1 In order to evaluate the percentage of PSAP 
staining over malignant prostate cytoplasm and the shape of 
the nuclei. PSAP protein, cytoplasm, and nuclei need to be 
labeled by three different dyes. In both cases, multiple dyes 
can be attached to several substances in a single structure or 
each dye can be attached to a specific structure. Although 



these dyes are visualized to have specific colors, they actu- 
ally have complex attenuation spectra. Therefore, in regions 
where the. stains overlap, the quantification of each stain 
component cannot be obtained ai any single wavelength be- 
cause the optical density ax any single wavelength or a broad 
range of wavelength is determined by the total attenuation 
from the multiple stains. 

Recent studies have shown that the optical mixtures of 
several absorbing substances can be analyzed in certain 
cases. Where the absorption bands of the individual compo- 
nent of the mixture do not overlap, the analysis is simple and 
straightforward. Wavelengths are selected where each com- 
ponent in turn displays an absorption adequate for measure- 
ment and the other components display negligible absorp- 
tion. Thus one attenuation measurement at each chosen 
wavelength suffices to determine each individual substance 
in the mixture. This technique has been successfully used in 
pathology specimens by Bacus era! 2 " 4 The advantages of 
ibis technique are that (I) it provides insight into the bio- 
chemical composition of the cell or tissue; (2) image acqui- 
sition is simply performed by using narrow band filters. Un- 
fortunately, current cytochemistry and histochemisuy ofTer a 
very limited number of such dyes with nonoverlapping spec- 
tra. 

During the past few years, the growing interest in color 
processing has led to (he development of spectral decompo- 
sition techniques using three channel color video cameras. 
The color video model uses three broadband filters in red, 
green, and blue to obtain a color image. Various color space 
transformation models from RGB have been developed. 5 "" 
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In the transformed three^imensional color space, color ob- 
jects can be distinguished. After a user de6nes o thresh Id, 
the image can be segmented The amount of each stain in 
objects with similar color can be determined by optica] den- 
sitometry. This technique allows object discrunijiation with- 
out prior determination of the spectral properties of the 
stains. Image acquisition is also Simply done using a com- 
mercially available RGB camera. The RGB technique is 
based on the human visual system with broad attenuation 
spectra. However, the spectral basis provided by the RGB 
detectors will in genera] be suboptimal for the specific stain 
combinations. Further, this technique docs not permit in- 
sights into the biochemical composition. The user defined 
thresholds may result Sn some errors in segmentation, which 
in turn affect the accuracy of quantification. Moreover, spec- 
tral information is lost during image acquisition by using 
broadband filters. Finally, spatial resolution is lower in RGB 
cameras than in monrxhroraic cameras because three pixels 
(red, green, and blue) must be interleaved, and the spacing 
between pixels is increased. 

To overcome the drawbacks of the above two techniques, 
we have developed a new color decomposition technique 
based upon the differences of mass absorptivity ratios at any 
two wavelengths. This technique can be used to quantify 
each stain component in a mixture or the stains in various 
structures. We Stan with a general formulation of the mul- 
tiple wavelength algorithm according to Laioben-Becr law. 
From this general formulation, dual and triple wavelength 
algorithms for two and three colon are derived. We describe 
the modifications that were necessary to implement this 
spectra! decomposition technique on a conventional trans- 
mitted light microscope. An excellent review of ths basic 
hardware and the hardware setup procedures for video mi* 
croscapy are provided by Incite** We also provide the de- 
tails of specimen preparation that are necessary for accurate 
measurements. Results of preliminary experiments are pro- 
sen ced 



IL THEORY 

A. Muttfpte wavelength algorithm 

In this discussion we assume a video microscopy system, 
where, to hist order, the pixel values of the digital image ate 
proportional lo the distribution of light intensity incident on 
the detector (/(xoO^Jr^x.y). In this system, when light of 
wavelength X ft passes through N components of concentra- 
tion cj(xoO and thickness r/x.y), we have from the 
Lambert-Beer law, 12,14 

C 4 (Ay) = C0 4 (x,v)exp| -£ m,^(x.y)/ 4 (x,y) j, (I) 

where m Jk is the mass absorptivity of the Jth component 
obtained at wavelength A*. C*(x,y) and C0 4 {x,v) are the 
images obtained through die N components and in their ab- 
sence, respectively. 



We define A k to be the absorbance at wavelength X A and 
define the mass density as ^(x,y)=c l (x i v) / i (x ? y). A t can 
then be written as 

If N sets of images G k (Jc= 1,2 N) are acquired at /V 

different wavelengths, X, ,Xa,...,X*,....X„, we obtain the fol- 
lowing N linear equations: 

A , (x . y ) = m , , $i (x.y) + m 2I Sfay ) + - - * + m«ri O0 * 

A 2 [x.y)-m (2J|(x v y)+m B 4(x.y) + ■ * ' +m>j**(x.y). 

Equation (3) can be written as a matrix equation: 



where we define 



(4) 













A* 






A = 










.A„. 








0I|| 


m 2 , 










■ i * 


• * . 






m iN mix 





m ffH t 

If \M | Q, we can solve for b\ as 
where 



(5) 





A, 


m 2 \ " 










At 




' Km 








• ■ • 


■ « * ■ ■ 




< 






A H 




m ftH 










A, - 


4 *w 










A 2 » 








A,= 






■ • •« 










A M " 


" m M< 






and 














m i\ 






-U 


a: 




m l2 






-u 


At 






* ■ • • 


mm m ■ 




mmm 




. m lff 


miff - 






Ah. 



(6) 



If the absolute value of the mass absorptivity, m l} {I J 
«1,2.3 N) is known for each of N components and at 
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each of the N wavelengths, then we sec from Eq, (5) thai the 
mass density distribution f each of the N components can 
be determined using the N absorbance images A k (* 
= 1,2,. JV) obtained at these same wavelengths. 

In practice, ii is difficult to determine the absolute mass 
absorptivity of the stain components in tissue. Although ac- 
curate measurements of each stain absorptivity spectrum 
could be made in solution, the stain absorption spectra is 
different in tissue than it is in solution. Thus, accurate mea- 
surements of the absorption spectra would be required in the 
tissue. In our current methods, we do not have an indepen- 
denl measurement of stain concentration in tissue. We have, 
therefore, chosen to measure the shape of the attenuation 
spectra for each stain in tissue and then solve only for the 
relative mass density distribution of each stain. To begin, we 
write Eq. (2) in the form 



Equation (4) is then written as the matrix equation 
A~RD\ 



a) 



(8) 
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m AW^V. 



(10) 



The matrix R contains the wavelength-dependent mass ab- 
sorptivity ratios T each component Equation (8) can be 
solved for the N mass densities as 



where we define 
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(12) 



Equation (1 1) demonstrates that the different stain compo- 
nents in the image can be separated based upon the mass 
absorptivity ratios of each component. Wheo the absolute 
mass absorpciviiies are unknown* we can select a wavelength 
for which all of the mass absorptivity coefficients are non- 
zero and use that as the Nth wavelength in Eq, (9). Then 
relative S values can be obtained according to Eq. (11). 
Therefore, a mixture containing N components can be de- 
composed oi N wavelengths with N(N- 1) mass absorptiv- 
ity ratios. 



B. Two-wnvelenfjth alportttim 

Consider a mixture containing two stain components, S { 
and At two chosen wavelengths, X c and X*. the 
wavelength-dependent mass absorptivity ratios of these two 
components are o=r n and /J=r M , respectively. The absor- 
bances at the two wavelengths are A and B respectively. 
Then we cao solve for £, and 6^ as 



4(*,y)< 



*(x,y)= 



A(x.y)-0-B[z,y) 

aBi]xj)-A(x t y) 
w 2 5-(c?-0) 



(13) 



C. Throfr-wavelength algorithm 

Similarly, consider an image that contains three compo- 
nents, 0^, and 4j. At three wavelength, \ a% X 6 , and X,., 
the wavelength*dependem mass absorptivity ratios of the 
first component are ori=r M and o 2 =r^ % those of the second 
component are ^[=r a , and A =r 22' and those of the third 
component are y=>r y . and Ji=r n respectively, The abscr- 
bances at the three wavelengths are >4, B. and C. respec- 
tively. The three components can be determined as 
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Optical Dentity 



(14) 

In practice, the doubly labeled tissue sections have three 
aitenufliion components; tissue itself and two stains. We as- 
sume chat the tissue section ilself appears translucent and has 
no color. We approximate che tissue attenuation ac any wave- 
length as A. In this case, we can replace the A and B in Eq. 
(U) by A r =.4-A andfl'=B-A. Similarly, the A, B, and 
Cin Eq. (14) can be replaced by A'- A- A and fl'-S-A, 
and C = C- A, respectively. 



III. METHODS 

A. Microscope hardware 

In order to implement the multiple wavelength decompo- 
sition algorithm, a microscope image analysis system was 
configured using; a transmitted light OLYMPUS BH-2 mi- 
croscope which is interfaced to a PC 486 DX3J MHz com- 
puter with 16 Mbyte of ram and $00 Mbyte disk and an SVG 
monitor. The necessary microscope functions, including 
translation stage positioning, light illuminfltion, and filter se- 
lection, are interfaced to the PC 486 computer using a MAC- 
2000 microscope controller and the PC RS232 serial pott 
The traditional microscope stage was replaced by a step mo- 
torized x/y stage. The resolution, reproducibility, and accu- 
racy of the stage are 0.1, 4 t and 3 ft respectively. Maximum 
stage movement speed is 30 000 /i/s. The focus was also 
controlled with a z axis stepper motor, with resolution of 
0.01 pl and movement speed of 27 000 ji& The filter wheel* 
with positions for six filters, was also motorized Tor com- 
puter control. The light intensity was either manually, or by 
computer control of a DC stabilized power source. 

An OLYMPUS achromatic objective* model Apo 20X, 
and a photo eyepiece with x5 magnification, giving a total 
magnification of xJOO. were used in image acquisition of all 
the image demonstrated in this paper. A XHXDC 1400 
monochromatic CCD camera (XTTJXX. Technologies Corp., 
Vancouver. British Columbia, Canada) (1280X1024X10 bit) 
was attached to the microscope photo tube and used as the 
detector. For oar system, with a total magnification of XJOO 
end the image matrix of 1280X1024. we obtained pixel di- 
mensions of 0.22 ix by 0.22 xl Image acquisition was under 
computer control. The sensitivity of the camera is approxi- 
mately 10"*- 100 lux (compared to the traditional 3-2000 
lux). The acquired Images were not Subjected to any pro- 



Fig. I, The line*/ correlation between meuurcd optica! dcnriiics md caU- 
bctttd sUMhnb. The coaetonoo ti io&ciied by Che slope 1.029. iotrrcrpi 
minai 0»009. and conelufon ceefficfem 0,999. 

ceiling algorithms other than the specific image decomposi- 
tion algorithms outlined in the theory section. 

Before the acquisition of any images, the light path was 
adjusted for Koehler illumination, as described by InoueJ 
The illumination was found to vary by about /% over the 
imaged field of view. Any deviation in illumination unifor- 
mity or variation in sensitivity as a function of wavelength is 
compensated by the normalization procedure of the measure- 
ment technique described in Sec. II Although the microscope 
MTF was nor measured, the performance of the entire system 
was judged to be quite adequate by an experienced patholo- 
gist. Tissue structures were observed with sharpness al- 
though the magnification and the resulting pixel dimensions 
are at the diffraction limit. For the study of this paper, the 
stains covered large regions, with moderately slow spatial 
variation and the measurements were primarily In the low 
spatial frequency portion of the system MTF. We refer the 
reader to the text by Inoue** for technical details pertaining 
to the hardware setup necessary for accurate quantitative 
measurements of stain intensity. 

The system linearity was measured using a stepped neu- 
tral density filler of known optical densities. The stepped 
filter had ten strips with gradually increasing density from 
zero to one and the densities were independently measured 
with an optical densitometer. For measurements, the zero 
optical density was used as the reference intensity. The re- 
sults are shown in Fig. 1. The relationship between the mea- 
sured and known optical density was OD tl ^ xmtJ "l.Q3 
OD^-0.009 with a correlation coefficient of r==0.999. 

Two interference filters centered at 480 (Xj and 570 am 
(X*) with 10 nra band widths were used to acquire the images 
for the dual wavelength algorithm. Three interference filters 
centered at 480 (X 0 )> 570 (X*), and 630 tun (X r ) with 10 nm 
bandwidths ware used to acquire the images for the triple 
wavelength algorithm. These filters are manually switched 
on the 61ter mount in the base of microscope. The wave- 
lengths were chosen to provide a difference in the mass ab- 
sorptivity ratios at any two wavelengths among stain compo- 
nents. At these wavel ngths, we found chromatic aberration 
to be insignificant. Airy image shifts due to filter misalign- 
ments were measured and corrected using an alignment grid. 
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A small wavelength-dependent focus adjustment was also 
necessary for the 630 nm filler but was not required for the 
other fillers. Image magnification was the same for all wave- 
' lengths used. 

B. Sampla preparation 

All the prostate tissue sections were cut with a thickness 
of 4 p. A two stain test sample was created by double stain- 
ing a prostate lissue section with PSAP-DAB (brown stain) 
and hematoxylin (blue stain). In this sample, the cytoplasm 
was sirongly stained by brown color and nuclei were 

. strongly stained by blue color. To test (he efficiency of the 
uipJe-wavclcngfh color segmentation, a second prostate tis- 
sue section was triply labeled by the PSAP-DAB (brown 
stain), hematoxylin (blue stain), and KeralinrPastRed (pink 
stain). In this case, the cytoplasm is partially stained with 
brown DAB and completely siaincd with pink while the nu- 
clei are stained with blue. Each stain will slightly precipitate 
in the connective tissue surrounding the prostate glands. This 
partial precipitation occurs throughout the sample and we 
refer to this weak stain in the connective tissue as the back- 

" ground. In order to estimate the wavelength-dependent ab- 
sorptivity ratio of each stain component and the attenuation 
of nonstaincd tissue, four sets of control tissue sections were 
stained by (1) PSAP-DAB (brown) only, (2) hematoxylin 
(blue) only. (3) Keraiin-FasiRed (pink) only, and (4) none of 
these stains. In each case, the control samples were pro- 
cessed in the same solutions as the multiple Stain test 
samples. A detailed description of sample preparation and 
control slide preparation is for a set of two-wavelength ex- 

. perimems has been published previously. 14,17 The prepara- 
tion technique for the three -wavelength experiments is a 
logical extension of these techniques. 

C. Estimation ol relative mos9 absorptMtiefi 

The tissue section background values are obtained from 
measurements on corresponding regions in the nonstained 
(issue section at each wavelength. After correction for this 
wavelength-dependent background, ihe relative mass absorp- 
tivities are estimated by making optica) density measure- 
ments in regions of Stain uptake on the singly stained slides. 
For the PAP-DAB and FastRed stains, regions are drawn that 
include only cytoplasm (no nuclei) within the gland area. For 
the Hematoxylin stain, regions are drawn to only include the 
stained nuclei. To obtain an estimate of the precision of the 
measurement, multiple regions arc measured. Wc report ihe 
standard deviation of the absorptivity taiios. 

Although not necessary for these experiments, the full 
spectra of the three stains used in this paper were measured 
using the singly labeled tissue sections from one control 
batch. The measurements were made using a ZEISS micro 
ph tometcr 400 and the resulrs arc shown in Fig. 2. Because 
of the spatial resolution of this system, the measurements did 
not exclude regions of nonstain upiake. and these curves are 
only approximations to actual stain spectra that would be 
required in this study. 




400 440 480 520 560 600 640 680 



WivdengtfKnm) 

Be. 2. A spectra] iramrnjawce plot of PAP-DAD. hematoxylin. and 
Kentin-FfcnRed slain prodoeu in prostate tiisoe sections. The data are ob- 
tained from ZEISS mioopftotomcter 400. 



D» Test of system accuracy 

Although it is difficult to obtain independent measure- 
ments of the amount of stain in a tissue section, we config- 
ured a set of experiments where ihe same set of regions on a 
slide were imaged during die various stages of staining by 
the multiple components. By carefully repositioning (he slide 
using our computer controlled xfy stage we were able to gel 
very accurately placed sets of multiple images at the same 
tissue locations. By measuring the optical density in essen- 
tially identical regions after each stain component was 
added, the differential optical density at each location was 
obtained. This differential optical density was used as the 
measurement standard for the corresponding stain. Ar the 
same time (in the same batch), singly labeled control slides 
were prepared for each stain used in the procedure. These . 
were also imaged ai the various preparation stages to ensure 
that the subsequent staining procedures did not modify slain 
density. The control slides were then used, as outlined above, 
to obtain the relative absorprjvities of each component at the 
imaged wavelengths and these were then used to segment die 
multiply labeled Slides, again as described above. The mea- 
surements obtained from die multiple wavelength segmenta- 
tion were then compared with those obtained from the direct 
measurements at the same locations. 

IV. RESULTS 

The assumption that the tissue section itself has no color 
has been validated in an earlier study, 16 where the absor- 



Table 1. Absorptiviiy raiios. duel wavelength. 
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bonce or a 4 p (hick section without any stains was found to 
be approximately 0.010. 0.011. aod 0.009 at 480, 570, and 
630 nm. respectively. Any errors to this estimate contribute 
only small errors in our measurements. Therefore, 
A '-A -0.010. £'=£-0.011, and C'=C- 0.009 are used 
in Eqs. (13) and ({4) instead of A, 0. and C. Because the 
absolute mass absorptivity is not known, relative 8 values 
were obtained by assuming that the m, d and m 2t> are equal to 
1 in Eq. (13) and m u and m u and m^. are equal to 1 in Eq. 
(14). 




Re. 3. (a) PftcKomjcrognph of the imaging field on Ac (result uuuc aide 
thai is tuned with PAP-DAB (brown) and Hematoxylin Ibtuel 
(nuanJficaiion=x J00V Image iovcsHy Is Inversely proportional to the &d- 
lorpdon of hemoioxylio and PAP stain, lb) The image of to) captwrd ai 
wavelength 450 ml (c) The image of (a> cupturctf ai 570 nm. Id) THe fceowq 
component 1100x4) image of {■), Image intensity a proportional to dw 
absorption of che PAP slain, te) Tbe bhie component 1 100*4) iraa*e ofiai. 



From the singly labeled tissue sections, v/c obtained an 
estimate of a equal to 1.71 for the brown only slide and p 
equal to 0.5] for the blue only slide (see Table 1). Images of 
a prostatic carcinoma slide stained with PaP-DAB (brown) 
and hematoxylin (blue) are shown in Fig. 3(a) was obtained 
with white light, linages captured ai wavelength 480 and 570 
nm are shown in Fig. 3(b) and Fig. 3(c). respectively. Ftom 
these images, we can see that the absorption spectra of the 
two components overlap at the wavelengths used. U is diffi- 
cult to separate them adequately in these monochromic im- 
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Tam-E 11. Ateorpliviiy ratio*, triple wavtlcngik 




Absorptivity mio SD 
#n,<480)/«i,[630) 


Absorpviviry isrio SD 


PAP-DAB 

HEM 

Ftaltcd 


a,=2.44, 0.049 
£,=0.67,0002 
y, a^.Qa^ 0.052 


a : =],42. 0.021 
ft- 1.39.01(02 
7y=4.& 0.09 1 



ages. After applying the dual wavelength method with the 
absorptivity ratios of <* equal to 1-71 and $ equal to 0.51 » we 
obtained the brown and blue slain component images shown 
in Figs. 3(d) and 3(e), The image intensities in Figs. 3(d) and 
3(c) are proportional to the mass per unit area of stain mate- 
rial on the slide so that image intensities of these figures are 

- inverted compared to the image intensities in Fig. 3(a). Fig. 
3(b). and Fig. 3(c). It is obvious that the blue stain is local- 
ized in the nuclei and the brown stain is in the cytoplasm and 
diey are completely separated 

From the singly labeled tissue sections, estimates of 
a, =2.44 and 1.42 for the brown stain and 0,=O.67 and 
ft= 1.37 for the blue stain and y, =4.03 and %=6.39 for the 
pink stain were obtained (sec Table UX Similarly, the images 
of the tissue stained with PSAP-DAB (brown), hematoxylin 
(blueX and Keraiin-FastRed (pink) arc displayed in Eg. 4. 
Figures 4(a)-4(d) art images of the same field captured at 
white light, 480, 570, and 630 nm, respectively. Again, it is 
difficult to separate the three components by using a conven- 
tional histogram threshold method because the absorption of 
these stain components overlap. Applying the triple wave- 
length algorithm by using the ratios listed above, we obtain 
the three component images as showo in Figs. 4(e)-4(g). In 
these images, the stains are segmented completely. The 
brown stain is only found in cytoplasmic area, the blue stain 

. is only on the nuclei, and the pin); stain is only in the cyto- 
plasmic area. 

The system accuracy results are plotted in Figs. S and 6 
for the dual and triple wavelength algorithms, respectively. 

V. DISCUSSION 

Wc have presented a multiple wavelength image decom- 
position algorithm, which can be used to obtain relative dis- 
tributions of multiple components with different attenuation 
spectra, This method allows the segmentation of tissue siruc- 
tures marked by different dyes or the decomposition of stain 
components in a mixture when these dyes have overlapping 
spectral absorption properties. Two wavelengths and two 
mass absorptivity ratios are required to separate two srain 
components while three wavelengths and six ratios arc re- 
quired to separate ihree slain components. In general. N mea- 
sures at N wavelengths and knowledge of N[N- 1 ) mass 
absorptivity ratios are required in the decomposition of an 
.image comoining N components. 

The mass absorptivity ratios at any two wavelengths can 
be simply estimated from attenuation measurements on sin- 
gly labeled tissue sections. To eliminate any errors that 
micbi occur from variation* in processing, these singly la- 



beled tissue sections should be prepared as controls and be 
processed in the same hatch in which the multiply labeled 
specimens are prepared. In this way, errors due to the stain- 
ing variations from one batch to another or from one tissue 
type to another can be greatly reduced. 

The multiple wavelength decomposition method should 
be useful in a wide variety of applications* where compo- 
nents with different attenuation spectra need to be separated 
for distribution determinations and quantitative analysis. The 
meihod overcomes the drawbacks of exclusive filtering 
techniques. 3-4 Compared with RGB color segmentation 
methods, our multiple wavelength method may facilitate the 
determination of insights into the biochemical composition 
in the tissue. Most important, it leads to the accurate quanti- 
tation of the relative distribution of cech stain component 
The use of narrow band filters limits the loss of spectral 
information in the image acquisition process. In theory, the 
use of the narrow band fillers and the use of optimally se- 
lected wavelengths should reduce the random error in the 
distribution measuremenu. 

Figures 5 and 6 give a good indication of the accuracy 
and reproducibility of the measurements made using this de- 
composition technique. Although die "true" measurements 
are made using the microscopy system, they are obtained 
directly for each of the sample points plotted. The estimated 
values are obtained using the absorptivity ratios measured on 
the control slides to decompose the values on the multiply 
stained slide. In each case the rest shows excellent correla- 
tion between the two types of measurements. These results 
are consistent with reproducibility measurements made pre- 
viously using the dual waveleogth method. 35 

Variations in section thickness as well as variations in the 
staining, tissue fixation, and processing can occur for differ- 
ent sample preparation batches. The effect of these variations 
can be significantly reduced by measuring the mass absorp- 
tivity using interna] controls for each of the study slides. For 
example, the variations of the color due to the physical and 
chemical conditions in the staining process can be eliminated 
by determining the relative wavelength absorption coeffi- 
cients using the singly labeled sections that are processed 
along with each staining procedure. Our preliminary study 
(unpublished observation) shows that die color of the stain 
components changes very little (coefficient of variation 
<0.051 from batch to batch. After extended periods of stor- 
age, a bleaching effect can occur, causing the tissue stains to 
fade. This effect was not observed over the short time period 
of this study. 

In another study, we have successfully applied this mul- 
tiple wavelength technique in the quantification of PSAP 
staining in doubly labeled prostate tissue sections. 27 In this 
study, die quantifications of PSAP (either intensity and ex- 
tern over prostate gland cytoplasm) were used to test the 
hypothesis thai the amount f PSAP in prostate cancer cor- 
relates to the patient survival. Normal glands that were posi- 
tively stained with PSAP were used as an interna) control for 
each of the multiply labeled tissue slides. The PSAP activity 
of the prostate tumor was evaluated by the ratio of PSAP 
intensity (or extend) In malignant glands versus PSaP tnten- 




01/16/2004 15:50,FM_613 230 6706 



OGILVY RENAULT 




(a) 




10 





FM. 4\ (a) Phocomiciograpf) of the imaging field on the prosuic oUwc slide 
thai "n soloed wirh PAP-DAB (brown ^ Hcmaioiylin (bine), and Kcnrirv- 
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tional ro the absuqHion of hematoxylin and PAP suin. ( b> The image or (a) 
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570 am. (d| The imags of (a) captures! at wauelength 630 nm. (c> The beown 
component (100x<ft im^e of (a). Image IntmsJty iff prwponional in the 
tbsoiplion of ifat PAP yam. (f) The blue component 1 100X0 unoge or fah 
(gj The pink component \\OOXJ) ima$e of (a). 
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FIG. 5, Accuracy of dual wavelenjtJi decomposition: $ystan neuurereenl 
(43 samples) versos sundanh. (a) The ttdmued amount of the PAP-DAB 
ruin- Here r«0.9W with slope 0.995 and wXctucp 0/MZ (b) U* eslmaied 
amount of hcBuuoiylm stain. Here n=0.998 wiih slope 0.97E and intercept 
0.013. 



sity (or extern) in normal glands presented in the same slide. 
We found significant correlation between the PSAP intensity 
and extent indices and patient survival. 27 We are continuing 
to study this correlation in a larger group of patients using 
the same technique. If our study jo the large population con- 
firms this hypothesis, it is likely that ihe multiple wavelength 
alg rithm may offer a new means for quantitatively evaluat- 
ing and specifying the treatment for patients with prostate 
cancer. 

The developed multiple wavelength algorithm is indeed 
worthwhile in quantitative histology and cytology since 
many prognose ically and diagnostically important substances 
in the tissue are commonly studied by labeling them with 
multiple stains. Numerous studies have addressed the quan- 
titative immunohislochemical analysis of proteins or other 
antigens.*-™' 25 We believe that this technique will be useful 
in the quantitative analysis of a variety of biological materi- 
als in multiply stained tissue sections for which die stain 
absorptivity spectra overlap. The technique permits decom- 
position and quantification of multiple stain components in a 
mature and opens a new way to perform functional studies 
in various fields where the Lambect-Beer absorption law 
h Ids.* The Lambert-Beer absorption law is generally valid 
except in some cases, such as when chemical interaction be- 
tween ihe stain components occurs. 1 - 1 * In this case. Hie pro- 
posed algorithm would need to be modified after the spectral 




Fee. 6. Accuracy or triple wavelength decamposfion. System measurements 
on 31 samples veraui suuxJartfe, fi) The esUmaicd amoum of PAP-DAB 
«ain. Here r^0S96 *i* slope 1.023 and misccpi -0.011. (b) The c»i- 
m>lcd arociw of ScmaiOxyUn sea*. Hew r«D.992 wilb slope J.039 and 
tmtnxpi -0.016. (c) The eauriated amount of Pasttad auin. Here r *0.972 
wilh slope 0.919 and intercept -0.002. 



properties of the post-interaction stain components are 
known. 

Image acquisition for the multiple wavelength technique 
is Jabor and time consuming, especially when the number of 
stain components is large. In order to solve this problem, an 
automatic microscopic imaging system has been established 
in our laboratory by utilizing a computer controlled stage, a 
filler wheel, and a power supply for controlling the illumi- 
nating light. By using this system, the user defined region of 
interest in a slide can be automatically scanned. The effi- 
ciency of the moltiple wavelength image acquisition has 
been impr ved by this updated system. 

Finally, this study has been performed without the ability 
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lo obtain th& absolute spectral absorption curve of ti\t stain 
components. With only a knowledge of ihc relative absorp. 
tion spectra, only relative measures were mode. Thus, the 
main goal of this paper has been to demonstrate <he multiple 
wavelength algorithm in multiple stain segmentation when 
only absorptivity ratios of the slain components are required. 
The wavelengths utilized in this study were chosen to give 
good separation of the stains which thai used in our experi- 
ment. After making accurate measurements of the complete 
attenuarion spectra it would be possible to select wave- 
lengths that wouW give better separation. Instead- we have 
limited our .study to wavelengths in a range thai seemed to 
give minimal image distortions due to wavelength-dependent 
imaging effects. 
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concentration and oxygen saturation in breast tissue 
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. Hear-mrxared (NCR) spectroscopic diffuse tomography has been used to soap the hemoglobin coma 
turn asd the hemoglobin cor/gen Gatumtion quantitatively in cuQuellke phantoms and *□ detei 
overage values U\ vwo. A series of phantom cslibratioas were performed to achieve quantita 
accurate images of the mbsorptiob end the reduced scattering coefficients at multiple optical wnvel ex 
A leaat-aquajces Gt was applied to ahsoxption-ooefficient images at multiple NIR wavelength* to c 
hemoglobin images of the concencradon and the hemoglobin orygen satw»tian- Objects of vo 
hemoglobin canamtnaon and ca^gen saturation within highly nm tf Ari og media were localize* 
imaged to within 2596 of actual values. Tha average hemoglobin concentration and oxygen 
ration of breast tissue was measured in vivo to two women voluntctaa. The potential application fi 
^fngnniM of breast bomcrfl b discussed. O 1959 Optical Society of America 
OCISaxies: 170,5270, llO.B080> W0-3S80, 170.6510, 170.69*0, 170-3B30. 



1. Introduction 

Recent studies have suggested that diffuse optical 
signals measured over the near-infrared (NIB) spec- 
tral range may provide sufficient information for the 
noninvasive mapping of hemoglobin concentration 
and oxygen saturation in biological tissues. 1 " 6 One 
driving force behind the development of NIB spec- 
troscopy and imaging has been the observed local 
increase in angiogeneeis, hence in blood volume, 
within cancerous tumors that corresponds to a local 
increase in the hemoglobin concentration, which has 
been estimated to exist at 4-times contrast with nor- 
mal tissue. 2 * 3 - 7 Breast malignancies have been 
shown to have a lower median oxygen pressure than 
the healthy breast A fi Some studies have pointed to 
a link between the oxygenation level in cancerous 
tumors and tumor characteristics. 10-12 For exam- 
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pie, Brizel et aL 10 found tumor oxygenation Ieve 
be related to the likelihood of metastases, wh« 
studies hy Okunieff et al. lx and Hoeckel et i 
showed a low-corygenation status of tumors to bi 
indicator of the effectiveness of radiation ther. 
In light of these findings it appears that noninva 
quantitative hemoglobin-concentration and oxyi 
saturation measurements have the potential to 
vide beneficial infbrmatiori for diagnostic 
therapeutic decisionmaking. Hence the realize 
of NIB spectroscopic diffuse optical tomograph} 
fers an opportunity to obtain this information 
only noninvasively but also spatially resolved an 
modest cost. 

In the diagnostic setting the ultimate goal of s 
an imaging modality would be to detect and diagi 
tissue pathologies an the basis of the spatially 
solved spectroscopic recovery of the optical abfr 
tion properties of tissue. These properties car 
mapped into physiological parameter distributi 
of interest, such as hemoglobin concentration 
oxygen saturation. Success in the spectxoso 
dimension been demonstrated in terms 
print measures of bulk-tissue response, leading 
single numerical indicators of functional tic 
characteristics. 1 ** Likewise, consider able effort 
been devoted to spatially resolving optical proper 
at a single wavelength by use of diffuse optical ixr 
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ing techniques- 13 -" However, ther have been, very 
few] attempts to combine spectroscopic data with an 
imaging method to produce spatially resolved images 
of the functional parameters that can be derived from 
measured spectra. 

Although optical images have been formed from 
diffusing light at multiple wavelengths in vivo * ex- 
isting methods are unable to recover quantitative im- 
ages of hemoglobin-related functional parameters. 
In breast tissue the primary absorbers are presumed 
to be oxygenated and deoxygenated hemoglobin, lip- 
ids, and water, 1 *- 18 whichhave significantly different 
absorption spectra in the near infrared. Quantita- 
tively mapping the absorption coefficient within 
breast tissue allows the contribution from each com- 
ponent to be approximated, and, in principle, images 
of the spatial distribution of all m^jor chromophobe 
concentrations can be obtained. This requires that 
data be recorded from a sufficient number of wave- 
lengths and processed by an appropriate image- 
reconstruction scheme to recover quantitative 
absorption-coefficient distributions. Previous re* 
search has suggested that quantitative optical- 
property reconstruction demands that the nonlinear 
relation between the optical-property distribution 
and the measured optical response be preserved. 1 *' 20 
To date, such techniques have not been applied to 
spectroscopic tomography data, and it remains an 
open question as to whether the accuracy required to 
derive quantitative images of hemoglobin concentra- 
tion and oxygen saturation can be maintained in the 
absorption-coefficient reconstructions at each wave- 
length. 

In this paper we combine the well-known concepts 
of optical spectroscopy with diffusa optical tomogra- 
phy to investigate the potential of spatially resolving 
functional hemoglobin parameters on the basis of 
noninvasrvely obtained, multisensor, inultiwave- 
length NIR measurements. Experiments were 
conducted on homogeneous and heterogeneous 
tissue-simulating phonfmng ofhuman blood mixed in 
a scattering medium that consists of Intralipid and 
water. Varying the blood volume and the oxygen 
pressure within these phantoms demonstrates the 
verification of the capability to quantify both bulk 
and spatially resolved hemoglobin-concentration and 
oxyg en-saturation levels. Experiments reporting in 
vivo bulk hemoglobin concentration and oxygen sat- 
uration are also described. The results suggest that 
spatially resolved quantitative maps of hemoglobin- 
related parameters should be possible with NIR dif- 
fuse optical tomography in the breast-imaging 
context This experience is the first report, to our 
knowledge, of quantitative NIR tomographic images 
of the hemoglobin concentration and the oxygen sat- 
uration obtained from multispectral experimental 
data with model-based image-reconstruction math" 
ods that preserve the important nonlinearities be- 
tween the functional property distribution and the 
measured optical response. 



Z Methods and Materials 

A- Data-Acqiisltlon System 
Our current data-acquisition system (Fig. 1) i: 
signed for cross-sectional imaging of breast til 
It consists of 16 source and 16 detector optical t 
positioned in a circular array. 21 (An in-deptt 
scription of an earlier version of the instrument* 
can be found in Ref. 21.) The optical fibers w 
the array are attached to a radial positioning syt 
that allows for variable-diameter fiber moves 
which provides direct-contact coupling with 
breast. The position of the circular array can all 
varied vertically to allow lor selective or xnul 
cross-sectional slices. 

A TiiBapphire laser OLexel Laser, Model 710) ' 
an argon-ion pump laser (Coherent, Model Innovi 
is used as a light source for our imaging ays 
The laser is tunable between 700 and 850 run by 
of a single optical setup. The laser beam is as 
tilde modulated with an electro-optic modti] 
(ConOptics, Model 360-80) that is driven by a 
MHz signal amplified by a 75-W rf amplifier (Ad 
fier Research, Model 75A25D). The laser bea3 
passed through a source optical fiber and dete 
sequentially through each of the 16 detector op 
fibers by means of multiplexing the fibers by use 
linear translation stago [Anorad Corp.). A ses 
translation stage is used to multiplex the laser b 
through each of the 16 source optical fibers, cres 
a total of 256 measurements. A phofcomultij 
tube (Hamamateu, Model R-928) is used for detec 
with a housing that contains a heterodyning cb 
(ISS, Model K209)> that allows for detection of 
signal amplitude and the phase shift at a lower 
quency of 1 kHz. 

B. (mage-Reconstruction Algorithm 

The 256 data points acquired for a single wavelez 
are reconstructed into an image of tissue optical p 
erties (absorption and reduced scattering coefficie 
by use of a finite- element-based calculation of 
diffusion equation. 22 The solution is obtar 
through a least-squares minimisation [minimize 
of ^ in Eq. (1), below] of the difference between 
measured fluence-rete 4>* intensity and phase at i 
detector i and the calculated fiuenoe-rate 4^ inter 
and phase determined from an estimated set of tir 
optical properties^ summed over the total numb* 
of observations: 

where the asterisk represents conjugation of 
complex-valued fluence rate. Newton's metho 
employed to solve the nonlinear minimization p 
lem iteratively; this problem results from the £i 
element solution of the diffusion equation for 
values of 4^ Further details an the reconstruc 
algorithm can be found elsewhere. 5 MW* 
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system, during patient studies* 



* C. Tissue Optical Properties 
The hemoglobin concentration can be mapped ap- 
proximately solely on the basis of the absorption co- 
efficient of tissue, as calculated at 800 nm, and an 
assumed level of water with other absorbers. He- 
moglobin oxygen-saturation mapping requires the 
calculation of the absorption coefficient at a mini- 
mum of two wavelengths. Although, approximate 
images can be obtained from two wavelengths, in 
many of our experiments wb chose to reconstruct 
quantitative images at three wavelengths or more to 
derive the hemoglobin-concentraticn and the axygen- 
saturation images. This method tends to minimize 
the effects of noise and can be used to account for 
absorption that is due to additional chromophores 
such as water and lipids. Additional wavelengths 
can be added, if necessary, to account far other ab- 
sorb era, such as proteins. 

We assumed that tie calculated absorption coeffi- 
cient Is a combination of absorption arising from wa- 
ter, oxygenated hemoglobin CHb-O*}, deoxygaaated 



hemoglobin (Hb-R), and lipids. In these exp 
ments, we neglected the absorption caused by oti 
tissue chromophores. For theoretical absorption 
efficients, we used the measured values from Wra; 
of. 35 for Hb~0 2 and Hb-R, from Hale and Quany 214 
water, and from Quaresima ti aI. M for lipids, G 
plots of the absorption coefficients for water, lipi 
Ht>0 2 , and Hb-R taken from the data provided 
these papers ere shown in Fig. 2. .For fhll-imj 
scans, we used three wavelengths: 760, B00 f t 
830 nm. In choosing these wavelengths, we selec 
750 nm to emphasize the large difference in abso 
turn coefficient between Hb-0 2 and Hb-R, 800 : 
because this is the isobestic point for Hb-0 2 c 
Hb-R, and 830 nm because Hb0 2 is more absorb: 
thfrn Hb-R at this wavelength- Thephotomultipli 
tube detector quantum efficiency decreases rapidly 
wavelengths longer than 850 nm, making detect 
at longer wavelengths impossible with the cum 
setup. 

Our imaging method and algorithm also pruvid 



01/16/2004 15:52 FAX 613 230 6706. 



OGILVY RENAULT 



H127 



i 




9 in r dians and the log of the intensity A o: 
measured transmitted light: 
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FSg. 2, Plat erf tho NIR absarptioQ spectra fas 1 water, lipids, Hb-O^ 
andKb^L Values far the absorptSea coefficient are displayed for 
anatomically relevant coaoftntratioos for breast tissue (31% water, 
57% lipid, lOjiM hemoglobin concentration at 60% oxj&m satura- 
tico)."^ 3 ^ Relevant values are estimated for Hb-O a and HthR 
from data bv Tmrobexj e* oX.? and valuer tor tbo water and the 
lipid content ere estimates from' eh e data gumnmiiiad by Woodard 
and White," 



quantitative mapping of the reduced scattering coef- 
ficient \l s * of the same region. 19 Measured quanti- 
tative values for the reduced scattering coefficient are 
compared with theoretical values by use of an equa- 
tion reported by Steveren at oJ. 2 * The scattering- 
coefficient images generated axe not used in the 
calculations of h^oglobin-concentration or hemoglo- 
bin oxygen-saturation level. 

D. Imagino Experiments 

Our initial experiments were focused on the verifica- 
tion of the capability to quantify correctly the absorp- 
tion coefficient for oxygenated and deosygenated 
blood (more specifically, Hb-O a and Hb-B) in a turbid 
medium. We deoxygenated the phantoms by flush- 
ing them with nitrogen gas. Verification of the de- 
oxygenation status was obtained by measurement of 
the partial pressure of oxygen (pO^) with a co- 
oximeter (Chiron Diagnostics, Model Rapid! ab 800). 
This device was also used to confirm the hemoglobin 
content of our blood samples in terms of the grams 
per liter of hemoglobin mass in blood volume; this 
value was used to convert from blood volume to he- 
moglobin concentration. We found the ccKsrimetzy 
readings to yield useful measurements of p0 2 as long 
as the b ample pressures were highftr than 20 mmHg, 
below which considerable variance appeared in tbe^ 
p0 2 as a result of tranqxutafcicn of the samples frxmr 
the phantom to the device* Nitrogen bubbling for SO 
xnin was apparently sufficient to deonygenate the he- 
moglobin effectively in the phantoms. Using our 
standard setup, we imaged homogeneous phantoms 
of water, Intralipid, and blood with 1 source and 16 
detectors. Applying the forward solution of the dif- 
fusion equation, we. itm^tjrely jBtte^ the compute-, 
taqps^o the. measured dateita ik^^ 
tis^n^oj^pal properties for a hirmo^nacnirmt^^ 
- Thife two-dimenaional fit was performed on the phase 



= 2 <&w - An - & w - aw 9 

where differences (shown in square brackets) 
used as an approximation of the first derivative : 
independent of the initial phase angle and so 
intensity, as described in Eq. (2). The term ; 
minimized on the basis of the comparative chanj 
the calculated difference in phase 6 and the log o. 
intensity A c between detectors i and i + I and 
difference In the measured values of 6° and A° for 
corresponding detectors summed over the comf 
set of M — 1 pairs of detector sites. 

Both lipids and water are significant NTH. abs 
ere 10 * 24 and affect the estimation of hemoglobin - 
centration. In addition, the spectra of 1 
compounds have small abeorbance peaks near 
nm, similar to Hb-R, which can lead to an unde 
tunation of the hemoglobin oxygen saturation, C 
Bequently, for phantom studies the ahsorbano 
water was subtracted from the calculated absorpt 
coefficient spectrum at each wavelength. For m i 
measurements a percentage of the lipid and the 
tar absorb ance spectrum that was based on the 
erage lipid and water content was subtracted fi 
the ahsorption-coefficient spectrum prior to perfb 
ing the least-squares fit to determine the hemoglc 
concentration and the hemoglobin oxygen saturat 
A goal of these experiments is to show that po 
by-point mapping within a tomographic cross sed 
of either breast tissue or tissuehke phantoms is : 
able. To image the hemoglobin concentration : 
the oxygen saturation in this manner, we applied 
finite-element method reconstruction algorithm 
data acquired at multiple wavelengths. For e 
wavelength the phase and the intensity data axe 
corded at 16 rirxurnferenrial detector locations 
each of the 16 source locations. An image of 
absorption coefficient and the reduced^cattering 
efficient is formed by use of the reconstruction al 
rithrxL Each different absorption image is tl 
considered as a mapping of tj^absorptian ooeffid 
for that wavelength. A least^qoar e g fiiis appliei 
^^pcwit in-theiniagB, r fee riahaltis .a pcfofcby^Jc 
gj^a^m^tip3Lpf the h^oj^obin^oncBntrati^ * 
Tt^erxi^ leveK' within breast tissue 

The absorption-coefficient and the redvu 
scattering-coefficient error bars displayed in Fx 
8-7 and 9 (below) were determined experiments 
by use of the standard deviation of repeated m 
surementa btained at each wavelength and bit 
concentration. The experimental uncertainties 
ported in Tables 1-8 for the measurement of 1 
molar absorption coefficient, the Hb-R and i 
Hb-0 2 concentrations, and the total hemoglol 
concentration and oxygen saturation were det 
mined by consideration of the worst-case aituati 
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TOO 7S0 6M 

Wavotangth M 

Pi & 9 Comparison of measured values and thoee predicted by 
Staveren H td* for tha reduced aeattering codMant m« a 
ph^tom containing O.Gtt Intralipid in water. Experimental es- 
timates oa-fl reached by means of fitting the data Gum a homofe- 
neous phantom to the fiaite-clamMit sn&del 

that results from applying the calculated experi- 
mental bit or for the absorption coefficient. 

3. Results 

A. Reduced Scattering Coefficient 

Although the reduced scattering-coefficient images 
are not used in the estimation of the hemoglobin- 
concentration or the hemoglobin oxygen-saturation 
level, tie ability to calculate these values quantita- 
tively is an important step in the validation of the 
imagine algorithm. The quantitative values for the 
reduced scattering* coefficient at several NTR -wave- 
lengths for a 0.5% Intralipid solation can be -verified 
by use of en equation reported by Staveren et al^ A 
comparison of the values calculated from the equa- 
tion of Staveren et el- and the measured data from 
our system is shown in Fig- 3. 

B. Homogeneous phantom Studies 

A preliminary verification, of the ability to quantify 
the Hb-0 2 and the Hb-R content in a tissualike phan- 
tom was accomplished through experimental deter- 
mination of the molar absorption coefficients of Hb-R 
and Hb-Oj at different wavelengths by the imaging of 
different concentrations of oxygenated and deoaygen- 
ated blood in a 0.5% Intralipid and water solution 
contained in a plastic-beaker phantom. Blood- 
volume concentrations were converted to hemoglobin 
concentrations by use of the hemoglobin molecular 
weight (64,500 g/mole) and the content values ob- 
tained from cOKarimetay measurements; the hemoglo- 
bin content of the blood was found to be roughly 156 
g/1 of hemoglobin mass in the blood volume. The 
experimentally determined molar absorption coeffi- 
cient 5a the slope of a linear regression that relates 
the absorption coefficient to the concentration. Hie 
results for Hb-R and Hb-O a at 750, 800, and 830 nm 
are shown in Pig. 4, with numerical values tabulat e d 
in Table 1- (The r2 value reported for the linear 
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Pig. 4. Plot of the experimentally measured absorption ooeffi 
versus the cod obd trntioD ofblood for three wavelengths: (a) I 
and (b) HW> 4 . Measurements wer* recorded from a hon 
neons scattering medium of 0,5% Intralipid in a 72-mm-diar 
cylindrical phantom- 



ragreasion ia the standard correlation caeffici* 
The oflfaet at 0% blood does not match the value* 
water absorption reported by Hale and Querry. 2, 

To investigate further the ability to distirigt 
Hb-0 2 and Hb-K within a turbid medium, we 
tempted to reproduce the infrared spectra of 1C 
Hb-0 2 and lOpM Hb-R (hemoglobin concentrat 
calculated from blood volume) mixed with 0.6% 
tralipid and water. The NIR absorption-coeffic 
Bpectra were again measured by use of data fro 
source and 16 detectors. The spectra are plotte 
pi^ 5 and compared with the expected values 
water plus blood. Again, the curves are relatr 
accurate except for an overall offset. 

The offset in the absorption coefficient appear 
be due to the translucent millimeter-thick pis 
walls of the beakers need In contain the tissue 
phantoms. To verify this hypothesis, we meast 
the spectrum of the water plus QJ6% Intralipid a 
tion in both a plastic beaker and a thin-walled bal] 
(see Fig. 6). The offset was not apparent in the 
constructed absorption-coefficient spectrum from 
thin-walled balloon^phantom data. 
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Tnbte 1. Comparison of fee E*p*r!m«mal Dots r miration of mo MoUr Absorption Coefficient (SIobo of a Unoor Rogrwton to ttio Data] wftb 

ValuasotWr^ora^andHaleandQueny* 



Molar Abaorptioo CoeffldftDt Offect (Value at Oonl Blood) 



Substaoce 


Wavelength 


VinytiaJ. 


Bxperunefit 


Hale and Querry 


Erperimaal 


J*»Va 


Hb-K 


750 cm 


0.00076 


040081 ± 0.00009 


0-0028 


0.0052 ±0.0001 


0.97 


802 run 


0,00046 


aoooso ±0.00011 


0.0020 


0.0057 ± 0.0003 


0.92 




833 cm 


0,00086 


o.oooa& ±o.oooob 


0.0035 


OX 048 2: 0.0001 


a9Bi 




750 cm 


0.OOD31 


O.00041 ± 0.00004 


O.D02S 


a0039± 0.0001 


0S3. 




802 am 


0.00046 


CL00042 ± 0X0004 


0.0020 


0.0044 ±0.0001 


0.94' 




833 cm 


0,00003 


aoocw? ± 0.00004 


0-0035 


0.0085 * 0.0001 


0.931 



Using the assumption that the observed offset is 
caused by the plastic wells, we applied a correction at 
each wavelength for a given plastic-beaker phantom- 
Then, using' the spectra for Hb-B, Hb-0 2 , and water, 
we performed a leasb-squarea fit on the corrected 
spectral data to approximate the concentration of 
each of the three components. For the data pre- 
sented in Pig. 5 the fitting results after correction are 
shown in Table 2. 




850 



ojnoo 



C. In Vm Homogeneous Measurements 
Similar testa can be performed with in vino meas 
ments of breast-tissue NIR spectra. Altitu 
breast tissue is heterogeneous, we' can realize 
proYimations to bulk values by assuming horn 
neous tissue and using the same measurer) 
scheme as for' the homogeneous phantom atu< 
An interesting comparison in terms of bulk ht 
globin concentration and oxygen saturation cai 
obtained by the imaging of the breasts of young 
older female subjects. By performing similar 1 
spectral scans of a 33-year-old woman and a 
year-old woman, we applied least-squares fibs to 
absorption-coefficient data to approximate the I 
hemoglobin-concentration and the average he 
globin oxygen-saturation levels in each case, 
measured spectra and the corresponding le 
squares curve fits generated from the calculi 
concentrations (numerical values are shown in 
hie 8) and the spectral values of Hb-R, Hb-Oj>, 
ter. and lipids are plotted in Fig. 7 for each autr 
The concentration levels of water and lipid v 
constrained to the average anatomical values 
31% water content and a 57% lipid content 
breast tissue, as summarized by Woodard 
"White. 17 These values for the average water 
lipid content are approximate and can vary sig 
icantly with age and between individuals. To 
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§ 0.0050 



| 0.0025 
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700 TfiO 000 BS0 
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<W 

Fig. 5, Companion of the absorption corfflripnt calculated from 
fo^sectedvaloOT with measured vahrafeblcwd^ and water 64 : (a) 
lOpX decnygvnnfced hemoglobin 2b a 0J6% latxmlipid sad water 
solntiaa and <b> lOpM oo^gecated fcezno^ltaina0£%Ic&Hllpid 
and water aolutton. 
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Fig. 6. Ocnipariaaa of the atooiption-ooemden t valued from I 
and Qoaay 1 * wttb the absorption ooefficaeat of water meaaun 
a pWtk beaker vunmj a tlnn-wallad hallocm. Not* that tha o 
ia cot apparent with the Qua-walled balloon, 
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TbMdX floruits of the Lea si-Square F«te for Cefculaiins the Hb-O a and &* Hb-R Content* of o Hotnooonooija Phantom of Known Comport 

IOilW HomogloNn Goncentratfon, 05% Intraflpltf, and Water 



Hb-O a 



Calculated 



Approximated Real 



Hb-O a coDcentxation (p-M) 8 .9 a: 1.6 

Hb-R conccBlratioa GiM) M ± 0.9 

Total Hb concexxtratjco ftiM) 10-6 = 0.9 

Oxygon saturation {%) 54 1 11 

Preset value of the water content 100% 



20.0 
0.0 

10.0 
100 
100% 



Calculated 



3.7 r 1.6 
8.9 A 1.0 
12.7 ± 0.8 
2B~16 
100% 



Hb-R 



Eange ofBea) 



0.0-2.5 
7.6-40.0 
10.0 
0-25 
100% 



count for this variation, we included in Table 3 the 
uncertainty introduced from a possible 50% varia- 
tion in the average water and lipid concentrations 
in addition to the experimental uncertainty, 

D„ Heterogeneous Phantom Studies 
To verify the ability to quantify correctly the ab- 
sorption coefficient of localized blood at three NTR 
wavelengths, we measured a series of different con- 
centrations of blood objects. A 72-mm phantom 
(plastic beaker) containing 0.5% Intralipid in water 
and 0.5% blood was used with a 25-mm object placed 
just off center. The concentration of blood in the 
object (a thin-walled deer plastic cylinder) was var- 
ied from 0% to 2.0%. The resulting images at 750, 
800, and 830 run are shown in Fig. 8, The 
absorption-coefficient value from the center region of 
the object is compared with the expected values cal- 
culated from the known hemoglobin (calculated from 
the blood volume by use of the value of 156 gr beano- 
giobin/l of blood) and water concentrations {based on 
values from Wrey ef gJ. 23 and Hale and Querzy 2 *). 
These data are plotted in Fig. 9, and the results of a 
linear regression on the data axe presented in Table 
4, . 

The 25-mm objects used in the heterogeneous 
phantom studies are thin-walled (200-p.m) clear plas- 
tic cylinders. The effect of the 200~pjn4hick cylin- 
der walls was evaluated through imaging 
experiments by use of identical contents for the object 
.and the background; these measurements ' produced 
no localized evidence of an absorption-coefficient in* 
crease. In contrast with the millimeter- thick walls 
of the exterior beaker, the 200- ^m walls of the cyKn- 



Tafito a Numerical Results of the U&st^Squarcs Fits far DeteonMno 
tho Hemoglobin Concentration and the Homoalcbb Oxygen Saturation 
It* Two Subjects* 





33-YearOld 


62-Yeaz^Old 




Woman 




Hb-Oa ameantzattai <u.MJ 


44 2 4 


15^:2 


Hb-R am centra ticn QjJd) 


1928 


7±1 


Total Hb cosoatitntroo (jiM) 


63 22 


22 32 


Qorgen aaturstbn {%) 


70 2 6" 


6727 


Preset values of tie wster am tent 


31* 


31% 


Preset values of the lipid content 


57* 


57* 



*Tue water and the lipid caamntratbna were constrained (a the 
approximate values displayed 



der appear to be negligible in comparison with 
experimental noise. 

Absorption-coefficient images at multiple we 
lengths can be used to generate an image of the 
moglobin concentration and the oxygen saturati 
The merit of this approach was tested with a serie 
Intralipid phantoms containing objects with vari. 
concentrations of oxygenated and deoxygena 
blood. Figure 10 shows a sample analysis of a 
suelike phantom. In this case en S4-mm-diame 
phantom consisted of a background of a fully oxyg 
ated lOjtM hemoglobin concentration in 0.5% 



aoso 
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Fig. 7* NIR spectra and graphical results of a Iaastrsqaarea fi 
Jete rCttin e the nemo^lobia anucenUuttoii and the b empale 
axygen-MWWdCft level for (a) a ^year-old woman and (b) a 
year-old woman. - 
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Fig. 6. Collection of abttuptios-coaffideot images lor a 25-mm einulu abject in a 72-mxn circular phantom. Hie background 1a 
blood in 0-6% Intralipid and water. The object In the first row naa no blood. Hue object in the second row contains 1% blood, that h 
third rov 1.5% blood, and that in the fourth row 29b blood. Column (a) is measured at 750 am, column <b) at 300 nxn, and oolnmn < 
830 nm. 



tralipid and water. His phantom contained two 
25 -mm circular objects embedded within the back- 
ground region. The object positioned off center to 
the right contained the exact same contents as the 
background, except that the object was bubbled with 
nitrogen gas to deoxygenate the blood. The second 
25-mm object, which contained a fully oxygenated 
40]lM hemoglobin concentration in 0-5% Intralipid 
and water, was positioned left of canter. Using a 



4L0150 




onooo I > i i - i .. » 

aooco 0X025 cnoso o«oo78 aoioo acis ojoiso 

LHootW Mum jitkm QMlTlrfaM fi {fttti*) 

Fig. 9. CamparJacn of the absorption coefficient taken team the 
center of the object ahown in Fig, 8 with the absorptinn-ooefficient 
values ealealaied from the expected values 0 - 14 for blood and water. 



least-squares fit, one obtains an image of both 
hemoglobin concentration and the hemoglobin ( 
gen saturation is formed from images at two wr 
lengths (750 and 802 mn), as illustrated in Fig. : 

4. Discussion 

The molar absorption-coefficient values calculi 
from the curves shown in Fig. 4 are similar to 
results reported byWray et ni. 23 for Hb-R and Hb- 
as shown in Table 1. Similarly, the results for 
absorption-coefficient spectra for Hb-R and Hb 
(Fig. E) are accurate whan plotted against establia 
values; however, an offset appears in both cai 
Further investigation led to the determination t 
this offset was due to the presence of the plastic t 
Iter used to hold the phantom solutions. Fhant 
experiments showed that the offset was not evid 



TmW* 4» UneaM^rrakm Remits from tho Data Plotted In Fig. 



Result 




Wavelength 




750 can 


802 am 


SSSnm 


Slope 


1J«0 


0.859 


0.670 


Oflaet 


OjOOO 


0.001 


0.000 


value 


0.972 


0-9M 


o^ai 



a A elope of 1 indicates a 1:1 correspondence between 
abgorptian^ngfTiriint values from the center of the 6Y$&± in 
images in Kg. and the abaarp^n-coefiBdenJt results calculi 
from tbe aapaaed values for hemoglobin and water.^ 
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Fig. .10. Abaorptioa-cooSdent, heionglohin-mncenfcrafaon, and hemoglobin oxygen^aafcurstum images of on otfect with an incrc 
hemoglobin concentration (off center to the left) find an object with deoaygenated Mood (off crate* to the right) in an axygeaated-t 
bac3cgrovmd: (a) Absorption, coeffident in inverse millimeters at 750 urn, (b) absorption coefficient In inverse millimeters at 802 no 
total hemoglobin concentration in micromolar, (d) Hb-O, in micron olar, fa) hemoglobin crygen saturation in percent. Tha plot* 
on»-dancnflional transects through (f) the image in (c), (tf) the image in (d) a (b) the image in (e). Ideal profiles are plotted for compar 



in a thin-walled balloon phantom (Fig. 6). This find- 
ing suggests that the offset ifi due to the effect of the 
beaker walla and that no such, effect should be ob- 
served for in vivo measurements. After we account 
fir the offset resulting from the presence of the bea- 
ter walls, the results match to within 20% of the 
known content (the oxygen saturation of the deoxy- 
genated blood wag known to be less than 26%) of the 
tissuelike phantom, as shown in Table 2. Ideally, 
future caKbnttioDs would be performed with thin* 
walled balloon phantoms; however, in practice. Oils 
makes it difficult to include heterogeneities and to 
bubble solutions with nitrogen gas. 

Breast-tissue absorption coefficient baa been 
shown to decrease with age. 28 The results, dis- 
played in Pig. 7, show a dramatic and expected 
change in the overall breast-tissue absorption coeffi- 
cient between the younger and the older female vol* 
unteers. Similarly, the calculated hemoglobin 
concentration is significantly higher in the younger 
woman- This result corresponds well with the 
known, on-average, decrease in breast-tissue density 
with age. 27 Quantitative results far the 62-year-old 
subject are similar to data reported by Tromberg et 
cd. 1 and Fantini et alfi The average hemoglobin 
oxygen-saturation values measured for the two 
healthy women of 68% and 70% axe nearly the same 
as the results for healthy breast tissue that were 
measured by Tromberg et at (66% and 68%). x 



The images shown in Fig. 8 highlight the abilit 
recover quantitatively the absorption coefficient 
olgeots within a scattering medium at multiple wf 
lengths. The absorption coefficient at the cente 
the object when compared with known values (plol 
in Fig. 9) shows a nearly one-to-one correlation. ' 
slope of a linear regression to the data illustrates t 
errors between 4% and 15% occur for the differ 
wavelengths examined. The quantitative assi 
ment of the absorption coefficient at multiple we 
lengths should lead to the ability to quantify 
hemoglobin concentration and the hemoglobin c 
gen saturation* 

For objects with known water content 
absorption-coefficient images at two wavdeng 
should be sufficient to image the hemoglobin cone 
trafion and the oxygen saturation. The hemoglol 
concentration and the hemoglobin oxygen-saturat 
images displayed in Fig. 10 correctly recover 
increased-liemoglobin^nc . object in 

hemoglobin-concentration image and the deoxyg 
a ted object in the hem globin oxygen-saturation 
age. The absarpticQa<oeffiaent images at 750 i 
800 urn are also displayed in Fig. 10; both, the 
creased hemoglobin concentration and the deoxyg 
ated objects are evident as Strang absorbers in 
750-nm image, whereas only the increaf 
hemoglobin-concentration object is recovered in 
800-nm image. (This result is expected — see 
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Hb-0 2 and the Hb-R spectral signatures in Kg. 20 
These results demonstrate an ability to quantify and 
to map the hemoglobin concentration and the hemo- 
globin oxygen saturation accurately, which may pro- 
vide useful diagnostic information. 

The results presented here have relied on the 
assumption that other absorbers play an insignifi- 
cant role in the absorption at NTR wavelengths. 
According to White et al., 1 * breast tissue can be 
assumed to be approximately GG% adipose tissue 
end 33% mammary gland for the average healthy 
adult human female. Using values from Woodard 
and White 17 demonstrates that the mammary gland 
has, on average, 51% water, 31% lipids, and 18% 
protein, whereas adipose tissue has an average 
composition of 21% water, 71% lipid, and 4% pro- 
tein (by mass). These values are approximate and 
vary significantly among individuals; however, on 
the basis of these data average breast tissue can be 
approximated roughly aa 31% water r 57% lipids, 
and S% proteins* 1817 Hence it appears that it will 
be necessary for in vivo studies to include calcula- 
tions of the lipid and the water contents and their 
contributions to NTR absorption. 13 Further study 
is required to determine the number of wavelengths 
necessary to extract accurately the absorption re- 
sulting from Hb-R and Hb-0 2 from the overall ab- 
sorption measured as a result of heme glob in, water, 
lipids, and other possible absorbers. 

In vivo estimates of the hemoglobin concentration 
and the oxygen saturation are currently being ob- 
tained with two to four wavelengths by several re- 
search groups that arc using NTR optical diffuse data, 
typically in the fbrm of bulk sample measurements. " 
Our research suggests that using Newton's iterative 
method of reconstruction baaed on a finite-element 
implementation of the diffusion-equation approxima- 
tion allows accurate hemoglobin-concentration and 
hemoglobin oxygen-saturation images to be obtained 
with KIR optical diffuse tomography, which can spa- 
tially discriminate localized changes in these quanti- 
ties. It is evident that bt vivo mapping of the 
hemoglobin concentration and the hemoglobin oxy- 
gen saturation by use of two wavelengths is quanti- 
tatively accurate only if the absorption caused by 
water and lipids is either insignificant or known and 
their content is constant throughout the tissue being 
examined . Heterogeneous m apping of the hemoglo- 
bin oxygen saturation by use of three wavelengths is 
likewise accurate only if the absorption cawed by 
water or lipids is either insignificant or a known con- 
stant. Neither of these assumptions appears to be 
accurate for breast tissue, which obEges the need for 
= znultispectral imaging to calculate concentrations of 
' all these relevant chromophorea simultaneously. 
Furthermore, using multiple wavelengths reduces 
the effects of noise and other errors when deriving 
these functional parameter maps from absorption- 
coefficient images, as was recently suggested by 
McBride etal.** 



5. Conclusions 

This research has confirmed the ability to a 
quantitative hemoglobin-concentration and 0x3 
saturation images (accurate to 15%) of embedde* 
jects within a, highly scattering media in w 
spatially resolved localized perturbations of t 
quantities are recovered with model-based in: 
reconstruction methods that preserve the nonli 
relations between the property map itself and 
measured multispectral response. In vivo tests 
demonstrated the potential to determine accur, 
the bulk hemoglobin concentration and hemogl 
oxygen saturation as well, with single-sot 
iUunination data sets. The effects of other abs 
era, such as local variations in water and 1 
content, warrant further investigation. Analys 
the composition of breast tissue, as taken from ] 
lished values, has suggested that these other ahs 
era will be important, which indicates the 
minimum of four wavelengths will be needed to 
dues quantitative images of tissue functional pax 
etere in the breast. This conclusion contrasts 1 
current practice in which typically data from only 
wavelengths are pursued. 

We are grateful to V. Quareairaa, S. J. Mate 
and M. Ferrari for providing the lipid-spectrum d 
We are also grateful to J. Prewitt for help with 
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Abstr«ct-A spectrophotometry mulUwavelength method for the simultaneous drtenniriallofl ^ ^ 
hSobln datives Including oxyhemoglobin, deoxyhemoglobin c^bo.yhcmo^ob.n f^^^^^ 1 " 
hTbren ^dcvdoocd The choice of analytical wavelengths was made by the method* of informauon factors 
Ed Et pro^Z^ZnJ^ foVX quantitative analysis wbs dcveJopeU ^^^X^^ 
mXrVi^^^amraio^ algebraic background correction, and combined melhoda Including linear pro- 
gramming and algebraic background correction. 
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Respiration, the source of cellular energy, is oxi- 
dation of organic substances by oxygen. Oxygen trans- 
port from lungs to tissues and reverse transport of 
carbon dioxide is performed by hemoglobin, the IVno- 
tional protein of erythrocytes. The function of the main 
hemoglobin derivatives including oxy- (HbC^ carboxy- 
(fTb(CO)A met- (metHb), and reduced hemoglobin 
(dcoxy hemoglobin, Hb) is determined by a ligand bound 
to the sixth coordination bond of iron atom of the 
heme moiety (O a . CO, H,0). Methemoglobin does not 
participate in the transport of gases but under certain 
pathological conditions, it can constitute up to 2(f/o of 
total hemoglobin. Investigation of the transport func- 
tions of erythrocytes requires the simultaneous moni- 
toring of the main hemoglobin derivadves. Ligand 
binding induces redistribution of electron density of the 
heme, influences magnetic moment and spin state, and 
causes optical activity. Electronic conformational in- 
teraction generates absorption bands In the visible part 
of the spectrum from 390 to 650 run that are charac- 
teristic for each derivative. This suggests that spectro- 
photometry can be used to quantify the main 
hemoglobin derivatives. 

Due to the complex nature of the subject, simple 
methods of sr^ectrophotometrical analysis lack required 
precision. According to ICushakovsky [1], hemoglobin 
concentrations arc assayed by transformed on into a 
system containing lesser number of components using 



* To whom correspondence should be addressed. 



chemical intcrconversion. Non-equivalence of the solu- 
tions before and after the reaction is a disadvantage of 
this method. Another solution of the problem is the 
method of Firordr. Disadvantages of this method in- 
clude complicated choice of analytic wavelengths due to 
similarity of absorption spectra of the assayed compo- 
nents, lack of individual absorption regions, and increas- 
ing error in case of systems containing more than four 
components. The problem can be solved by redetermined 
systems of linear algebraic equations (SLAE). 

The goal of the present work was CO develop a 
spectrophotometric method of simultaneous assay of 
the main hemoglobin derivatives including oxy-. deoxy-. 
carboxy-, and methemoglobin* 



MATERIALS AND METHODS 

Samples of heparin-treatcd Wood were withdrawn 
from the elbow vein of donors. Blood solution (1%) 
was prepared by mixing 20 ml of distilled water. 6.06 ml 
of 0.04% ammonium solution (lo clear the solution), 
and 0.5 ml of the blood; 1-2 min later (after the 
hemolysis), 25 ml of 0.0667 M potassium-sodium 
phosphate buffer (pH 7.2) were added. Total volume 
was adjusted to 50 ml with distilled water. Spectra were 
recorded at 450-650 nrn in a quartz cell (1 can thick) 
using a SF-20M double-beam spectrophotometer with 
automatic aperture regulation. Assays were performed 
within 1 b after blood withdrawal. 
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RESULTS AND DISCUSSION 

Quantitative spec troph otometrlc analysis is reduced 
d wn to solution of the system of Firordt equations: 



0) 



where i = 1, 2, m; m is the number of analytical 
wavelengths; n corresponds to the number of quanti- 
fied components; e) ! are the molar extinction coeffi- 
cients of component j at wavelength Cj is the quan- 
tified concentration of component/ corresponds to 
absorption at wavelength V 

The methods of solution of redetermined SLAE osed 
in the present work are considered below. 

The least aqnare method (LSM). This method is the 
classical solution of redetermined SLAE. The criterion 
of the solution u the minimum of the sum of squares 
of deviations of the experimental spectrum from the 
expected one: 



V - 



f « V 
A** - X A c. 



= mm. 



Solution of system (I) by the least square method 
can include negative concentrations of certain compo- 
nents that lack physical sense. Negative concentrations 
are usually associated with similarity of absorption 
spectra of certain components. Moreover, negative 
concentrations appear when the matrix contains many 
components some of which are present at low concen- 
trations or are absent. Negative values m solution of 
the system (1) can be excluded in two ways. 

I. The following sum is used as the criterion of 
optimal solution: 



P) 



where cy arc the negative elements of the concentration 
column vector. 

2. The values can be obtained by iteration so that 
the negative concentration* obtained during the first 
iteration are assumed to be zero and the column cor- 
responding to these components are excluded from the 
matrix. 

Algebraic background correction method (ABC). In 
systems with unknown qualitative composition, absorp- 
tion is expressed as follows: 

/' = Z£|ic J + I fl X (4) 
BIOCHEMISTRY (Moscow) VoL 63 No. 10 1998 



This redetermined SLAE is solved by the least 
square method so that the first n columns of the matrix 
of the system coefficients c ntain the e] ( values and 
the subsequent / + 1 columns contain the X, values. 
To lower the calculation error due to rounding, the 
X values in Eq. (4) are replaced with the X - X values 
(where X is the mean wavelength of the interval) and 
for equidistant wavelengths, with their consecutive 
numbers. The most important factor determining the 
validity of this method is the difference between the 
absorption spectra of contaminating and main com- 
pounds. Efficiency of the ABC method increases when 
the impact of the main compound increases (polyno- 
mial order k > I + 1 . 

Linear programming method (LPM). The LPM is 
used to exclude the effect of background absorption 
during quantitative analysis of a muhicomponent mix- 
ture, in this case, absorption equals: 



A x ' * = Z $ cj 



(5) 



The system requires the minimal differential spectrum 



mm 



under constant £ x A 1 this is equivalent to the condition 
max ^ j^ef c r This means that the maximum of the 
function 

Z = X Z c> q 

should be found constraining A** - jE r)' cj > 0. 

According to the LPM. the solution of the quan- 
titative analysis problem corresponds to that of the 
method of analytical bands selected in the most clear 
region of the spectrum of contaminations but in the 
case of LPM, the spectrum of the contaminations d es 
not have to be known [2]. Differential spectrum cor- 
responds to the spectra of non-Identified components. 
The simplex optimization method is used to solve system 
(I) that includes limitation of negative concentrations. 

Combined methods including LPM and ABC. If tho 
spectrum of the mixture lacks individual absorption 
regions of the main compounds this can be used to 
analyze wavelengths corresponding to minimal absorp- 
tion of contamination. Such wavelengths can be se- 
lected by the LPM. In another words, the LPM selects 
analytical wavelengths with minimal or zero absorp- 
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lion of ooatamination. Tbe LPM solution is correct 
only if the studied spectral region includes wavelengths 
with zero background absorption. Otherwise, the 
method does not eliminate the influence of contami- 
nation on the results but only minimizes such influ- 
ence. 

Thus, it is of interest to combine the LPM and 
ABC algorithms in a single method. The first combi- 
nation is the selection of analytical wavelengths with 
minimal absorption of contamination by the LPM and 
subsequent polynomial approximation of tho back- 
ground spectrum at these wavelengths by ABC. 

The second possibility is to set np a problem of 
linear programming possessing ABC properties. The 
criterion of the solution of this problem is: 



U = Z - £ e> c. X flAi = 

/■J y J°t ' 3 J 



mm, 



(6) 



and constraints include the following expressions: 



I ati > 0; Cj > 0; - X £> c, - X aX > 0. 



LPM + ABC combination enables analysis of a 
mixture by approximating the background absorption 
with polynomials of higher exponent versus the spec- 
trum of tho main compound. 

The developed programs were tested for mixtures 
with known concentrations of the components. Solu- 
tions of methyl orange (existing in two forms: acidic 
and alkaline) were used as the reference. The data 
obtained by various methods are compared in the table. 
The results are in good agreement with the reference 
values, thus suggesting efficiency of the programs and 
their possible application in a more complicated sys- 
tem, i.e., hemoglobin derivatives. 

Gas composition of the whole blood solution was 
studied by computer programs used for analysis of 



t ■ 10~ 3 , HUr/mola per cm 
18 




Fig. 1. Absorption spectra: J) oxyhemoglobin; 2) dcoxy- 
hemosJohfa; J) aufcojyHoiictalobfci; 4i methonogjottn [3J. 



absorption spectra of the multicomponent systems of 
hemoglobin derivatives containing HbO>, Hb, HbtCO)*, 
metUb. Molar extinction coefficients used for analysis 
were taken from [3] (Fig. 1). Analytical wavelengths 
were selected by the information coefficient method f4J. 
To evaluate the optimal selection of wavelengths and 
error of analysis, the information coefficient matrix is 
constructed from the following elements: 



Z et 1 



(7) 



The Ry element indicates the fraction of informa- 
tion of tho / band on tbe determination of the J com- 
ponent of the mixture. The higher are the information 
coefficients, the lower are the errors of the solution. 
Also, to determined analytical wavelengths the differ- 
ences - ef* and ratios e> / ej* were used to select 
the points of maximal differences in oxyhemoglobin 



Comparison of alkaline and acidic fractions (%) in methyi orange mixtures by LSM, ABC, JLPM, and LPM 
ABC 



Mixture 
number 


Reference 


LSM 


ABC 


LPM 


LPM + ABC 


alkaline 


acidic 


alkaline 


acidic 


alkaline 


[ addle 


alkaline 


acidic 


alkaline 


acidic 


1 


6 


94 


5.9 


94,1 


5.8 


9da 


5.8 


94,2 


6.1 


93.9 


2 


9.7 


90,3 


9.2 


90.8 


8.9 


91.1 


9.3 


90.7 


9.5 


90.5 


3 


88.6 


11.4 


88.5 


11.5 


88.8 


11.2 


88.5 


IU 


88.6 


U.4 


4 


94.1 


5.9 


94 


6 


94.6 


5.4 


93.8 


6.2 


94.3 


5.7 


5 


96.6 


3.4 


96.4 


3.6 


96.6 


3.4 


96.5 


3.5 


96JS 


3.4 



1 
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Fit* 3L Eflcd of tnt«i hemoglobin concent ration dctcrmin**! 
by the rtandtrd cyAitmcthcntoglobin melbod on the sum or 
cononntritlofa of hemoglobin derivative* calculated by ibe 
developed method. 



and carboxyhemogJobLD spectra (their spectra are quite 
similar thus contributing significant error in determi- 
nation of car b oxyhemoglobin concentration). Then, 
these wavelengths were studied by LPNL Their com- 
position varied during analysis. Those wavelengths that 
had no effect on the mult were excluded from the 
experiment Minimal difference between experimental 
and reduced spectra at a certain wavelength wag an- 
other criterion of selection. 

To illustrate the developed method, the effect of 
nicotine intoxication in men on concentrations of 
hemoglobin derivatives was studied. The data of Fig. 
2 indicate that the body undergoes negative changes 
including lowering of oxyhemoglobin content due to 
increase in carboxyfaemoglobin concentration and hence 
lowering of hemoglobin oxygenation. 

The precmon of the method of quantitative analy- 
sis of mixtures of hemoglobin derivatives was evalu- 
ated. The sum of concentrations of hemoglobin deriva- 
tives calculated using the method would be equal to 
total hemoglobin concentration in the same blood 
sample determined by the standard cyanmethexnoglobin 
method; this was the criterion of correctness of calcu- 
lations of concentrations of hemoglobin derivatives. 
Effect of total hemoglobin concentration on the calcu- 
lated sum of concentrations of hemoglobin derivatives 



is shown in Fig. 3. Slope of the approximating line is 
0,997 and significance of appicodmation is 0.987, thus 
confirming the identity of the results. 

Thus, the multiwavelength method of determina- 
tion of the main hemoglobin derivatives ia developed 
and the software is created. The method can bo used 
in quantitative spectrophotometric analysis of any 
chemical systems with overlapping absorption spectra 
of the' components. A special multiwavelength module 
with fixed selection of analytical wavelengths can be 
constructed on the basis of the method and used in 
clinical practice. 

The work was supported by research programs 
''Universities of Russia." and "Siberia". 
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